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ABSTRACT 



> 

■ We present high spatial resolution integral-field spectroscopy of 28 elliptical (E) and lenticular 

\^ I (SO) galaxies from the SAURON representative survey obtained with the OASIS spectrograph 

. during its operation at the CFHT. These seeing-limited observations explore the central 

8"x 10" (typically one kiloparsec diameter) regions of these galaxies using a spatial sampling 
four times higher than SAURON (0'.'27 vs. 0'.'94 spatial elements), resulting in almost a factor 
of two improvement in the median PSF. These data allow accurate study of the central regions 
O ' to complement the large-scale view provided by SAURON. Here we present the stellar and 

gas kinematics, stellar absorption-line strengths and nebular emission-line strengths for this 
c/3 ■ sample. We also characterise the stellar velocity maps using the 'kinemetry' technique, and 

derive maps of the luminosity-weighted stellar age, metallicity and abundance ratio via stellar 
population models. We give a brief review of the structures found in our maps, linking also 
to larger-scale structures measured with SAURON. We present two previously unreported 
, _ kinematically-decoupled components (KDCs) in the centres of NGC 3032 and NGC 4382. We 

^ . compare the intrinsic size and luminosity-weighted stellar age of all the visible KDCs in the 

full SAURON sample, and find two types of components: kiloparsec-scale KDCs, which are 
older than 8 Gyr, and are found in galaxies with little net rotation; and compact KDCs, which 
have intrinsic diameters of less than a few hundred parsec, show a range of stellar ages from 
0.5 - 15 Gyr (with 5/6 younger than 5 Gyr), are found exclusively in fast-rotating galaxies, 
and are close to counter-rotating around the same axis as their host. Of the 7 galaxies in the 
SAURON sample with integrated luminosity-weighted ages less than 5 Gyr, 5 show such com- 
pact KDCs, suggesting a link between counter-rotation and recent star-formation. We show 
that this may be due to a combination of small sample size at young ages, and an observational 
bias, since young KDCs are easier to detect than their older and/or co-rotating counterparts. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: formation 
- galaxies: kinematics and dynamics - galaxies: structure - galaxies: ISM 
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1 INTRODUCTION 

In the era of the Hubble Space Telescope (HST), early-type galaxy 
nuclei (central 10-50 pc) have become a key component in the mod- 
ern picture of galaxy evolution. Empirical relationships between 
the central luminosity profile and numerous large-scale galaxy 
properties, such as isophotal shape, rotational support, and stel- 
lar populations, strongly suggest that the formation and evolu- 
tion of the central regi ons is in extricably linked to that of the 
galaxy as a whole (e.g. iFaber et a l. 1997). The ubiquitous pres- 
ence of supe r-massive black holes (SMBHs) at the centres of 
galaxies (e.g. iFerrarese & MerritJ bOOOt iGebhardt et alJIlOOd) is 
thought to be an important aspect of this relationship, control- 
ling global star-formation and dictating the central density distri- 
bution (e.g. Silk & Rees 1998: Granato et aL.2004: .Binnev..2004« : 
iDi Matteo. Springel & Hemauisj2005i) . 

Most observational studies of early-type galaxy nuclei have 
so far been based only on HST imaging, with a small number of 
objects also observed spectroscopically from space, either with a 
single aperture or a long slit. The limitations of these data be- 
come apparent when considering the complexity of structures found 
in the centres of these otherwise 'boring' objects. Spectroscopy 
is required to establish dynamical black hole mass estimates, and 
it has been shown that a single long slit i s grossly insufficient 
to constrain fully general dynamical models iVerolme et aljEoo3 : 
ICapDellari & McDermid 2005). To understand the stellar orbital 
structure of galaxy nuclei therefore requires information across a 
sizeable two-dimensional field, obtained with the highest possible 
spatial resolution. Additionally, HST imaging reveals highly com- 
plex gas and dust distributions (e.g. 'van Dokkum & Franxlll995l : 
ffomita et al. 2000; Lauer et al. 2005), which are sensitive indica- 
tors of dynamical torques that are not detectable from the stellar 
distribution alone. Understanding the kinematic and chemical prop- 
erties of such irregular gas distributions requires two-dimensional 
spectroscopy. //5r imaging studies also find broad-band colour gra- 
dients, indicative of distinct star-formation histories towards the 
very central regions (e.g .lCarollo & Danziget 1994b; Carollo et alJ 
ll997atlLauer et all2005l ; lcote et aU2006,) . However, to accurately 
quantify the populations requires spectroscopy on comparable spa- 
tial scales, of which there are rather few studies of early-type 
galaxy samples, due to the high signal-to-n oise ratio required (e.g. 
lKrainovic& Mf3l2004 ISarzi et aljEoo3) . Moreover, these and 
comparable g round-based studies looking at central properties (e.g. 
ICarolIo & D anziger 1994b; Trager et al. 2000; Thomas et al. 2005) 
are restricted to either aperture measurements or long-slit profiles, 
making the association between population changes and other struc- 
tural variations difficult to determine. 

In response to these issues, we have initiated a program of high 
spatial resolution integral-field spectroscopy of early-type galaxy 
nuclei. We select our sample from galaxies observed as part of 
the SAURON survey: a study of 72 representative nearby early-type 
galaxies and spiral bulges observed with SAURON, our custom-built 
panoramic integral-field spectrograph mounted at the William Her- 
schel Telescope (WHT) on La Palma (.Bacon et al.„ 2001. Paper I). 
The aims of this survey and the observed galaxy sample are de- 
scribed in'de Ze euw et alJ i2002l Paper II). For the sub-sample of 
48 E/S Os we present the stellar kinematic maps in Emsellem et al. 
{20041 Pap er III); the two-dimensional ionised gas properties in 
ISarzi et sl\ 12006. Pap er V); and ma ps of the stellar absorption line 
strengths in! Kuntschn er et all i2006l Paper VI). We refer the reader 
to these papers (and references therein) for discussion of the global 
properties of our sample galaxies. 



To maximise the field of view, the spatial sampling of SAURON 
was set to 0'.' 94x0'.' 94 (per lenslet) for the survey, therefore of- 
ten undersampling the typical seeing at La Palma (0'.'8). This does 
not affect the results for the main body of early-type galaxies, as 
they generally exhibit smoothly-varying structure on scales larger 
than the seeing. However, the central regions are poorly resolved in 
the SAURON data, preventing accurate measurements of the central 
properties. 

For this reason we have observed the centres of a subset of 
the SAURON sample using the OASIS integral-field spectrograph, 
mounted on the 3.6-m Canada-France-Hawaii Telescope (CFHT), 
Hawaii. This instrument allows optimal sampling of the natural see- 
ing on Mauna Kea, resulting in a significant improvement in spatial 
resolution. To ensure a high degree of complementarity between 
the OASIS data presented here and our SAURON observations, we 
configured OAS I S to have an almost identical wavelength domain 
and spectral resolution. The OASIS data can therefore be easily 
combined with SAURON data, e.g. for constraining dynamical 
models (e.g. Shaoiro et alj|200'^ . All data presented in this paper 
will be made available through the SAURON web page.^ 

The OASIS observations themselves reveal a number of pre- 
viously unresolved features in the centres of our sample galaxies. In 
particular, we find several galaxies which exhibit very young cen- 
tral stellar populations that can be directly associated with peculiar 
(misaligned) components in the stellar rotation fields. We explore 
this further by considering all objects in the SAURON sample with 
clear kinematically decoupled components, and find a link between 
the size of these components and their stellar populations. 

The layout of the paper is as follows: Section |2| describes 
the instrument, the observations and the basic data reduction; Sec- 
tion |3| describes the application of our analysis tools used to de- 
rive the stellar kinematics, gas properties and line-strength mea- 
surements (including the application of stellar population mod- 
els), and presents detailed comparisons with equivalent quantities 
from SAURON; Section|4|presents the two-dimensional maps of all 
measured quantities, and discusses their main features; Section |5| 
presents a discussion of young stellar populations in our sample, 
and the connection with KDCs; and Section|6|concludes. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 The OASIS Spectrograph 

OAS I S is an integr al-field spectrograph based on the TIGER con- 
cept {Bacon et aljl 995) and is designed for high spatial resolution 
observations. It is a multi-mode instrument, with both imaging and 
spectroscopic capabilities, and can be assisted by an adaptive op- 
tics (AO) system. OASIS operated at the Cassegrain focus of the 
CFHT between 1997 and 2002 (optionally with the PUEO AO sys- 
tem), and was transferred to the Nasmyth focus of the WHT behind 
the NAOMI AO system in March 2003. All observations presented 
here were obtained at the CFHT, hence in the following we solely 
refer to the OASIS CFHT configurations. 

Since most of the objects in the SAURON sample do not have 
a bright, nearby guiding source (mv < 16 within ~ 30"), the //8 
(non-AO) mode of OASIS was used. A spatial scale of 0'.' 27x0'.' 27 
per lenslet was selected for the observations to properly sample the 
generally excellent seeing at Mauna Kea, providing 1009 individual 
spectra in a 10" x 8" field-of-view. OASIS was configured to give 

^ http://www.strw.Ieidenuniv.nl/sauron 
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Table 2. Summai'y of observations and PSF estimates. 



Figure 1. SAURON survey sample of ellipticals and lenticulars. Objects are 
separated as 'field' (circular symbols) and 'cluster' (triangular symbols) as 
defined in Paper II. Filled symbols indicate the subset of objects observed 
with OAS I S at the CFHT. 



Table 1. Summary of the three OAS I S observing runs. 



Run 


Dates 


Clear 


1 


24-27/03/2001 


3/3 


2 


8-13/01/2002 


2.5/3.5 


3 


6-10/04/2002 


4.5/5 



Note: The dates indicate the duration of the run, although some of the nights 
were shared with other programs. The last column indicates the effective 
number of clear nights from the total number of nights dedicated to this 
program. 



NGC 


Type 


Run 


Tea; 


Source 


Seeing 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


1023 


SBO-(rs) 


2 


2 X 3600 


F555'W 


1.41 


2549 


SOO(r)sp 


1 


1800, 3600 


F702'W 


0.91 


2695 


SABOO(s) 


2 


2 X 3600 


Acq. 


0.91 


2699 


E: 


2 


2 X 3600 


F702'W 


1.67 


2768 


E6: 


3 


2 X 3600 


F555'W 


0.89 


2974 


E4 


2,3 


2 X 3600 


F555'W 


0.99 


3032 


SABOO(r) 


2 


2 X 3600 


F606'W 


0.92 


3379 


El 


1 


2 X 3600 


F555'W 


0.91 


3384 


SBO-(s): 


1 


1 X 2700 


F555'W 


0.72 


3414 


SO pec 


3 


2 X 3600 


F814'W 


0.86 


3489 


SABO+(rs) 


3 


2 X 3600 


F814'W 


0.69 


3608 


E2 


1 


1800, 2 X 2700 


F555'W 


1.13 


4150 


SOO(r)? 


2 


3 X 3600 


F814'W 


2.15 


4262 


SBO-(s) 


3 


2 X 3600 


F475'W 


0.60 


4382 


SO+(s)pec 


2 


3 X 3600 


F555'W 


1.11 


4459 


SO+(r) 


3 


2 X 3600 


F814'W 


1.53 


4473 


E5 


3 


2 X 3600 


F555'W 


0.80 


4486 


EO-l+pec 


1 


3 X 2700 


F555'W 


0.57 


4526 


SABOO(s) 


1 


1800, 3600 


F555'W 


0.94 


4552 


EO-1 


1 


2 X 2700 


F555'W 


0.67 




tj 
JJ, 


-2 
J 


Z X oOUU 


/ uz W 


U. /U 


4621 


E5 


3 


2 X 3600 


F555'W 


0.86 


5198 


El-2: 


3 


2 X 3600 


F702'W 


0.84 


5813 


El-2 


1 


3 X 2700 


F555'W 


0.87 


5831 


E3 


3 


3600, 1200 


F702'W 


0.95 


5845 


E: 


3 


2 X 3600 


F555'W 


0.91 


5846 


EO-1 


3 


2 X 3600 


F555'W 


0.58 


5982 


E3 


3 


2 X 3600 


F555'W 


0.77 



similar spectral coverage and resolution as SAURON, resulting in a 
wavelength range of 4760-5558 A, with a resolution of 5.4 A full- 
width half-maximum (FWHM), sampled at 1.95 A pix~^. Given 
the generally large velocity broadening in the central regions of 
early-type galaxies, this configuration is suitable for measuring 
stellar kinematics in these objects, and also covers key absorp- 
tion/emission features such as Mg b, a number of Fe lines, H/9, 
[OIII]AA4959,5007, and [NI]AA5 198,5200. 



2.2 Observed Sample 

In the chosen configuration, OAS I S has a spectral resolution around 
20% larger than that of SAURON, effectively increasing the min- 
imum velocity dispersion that can be reliably measured. For this 
reason, we restricted our sample to only the E/SO galaxies from the 
SAURON survey (which tend to have higher dispersions than the Sa 
bulges), and gave priority to objects with a central velocity disper- 
sion (as measured from Paper III) larger than the instrumental res- 
olution of OASIS (croAS ~ 135 km s~^). We observed a total of 
28 galaxies from the parent sample of 48 during three runs between 
March 2001 and April 2002 (TableQ. Table|2llists the galaxies and 
gives a summary of the observations. Fig. Q shows how these ob- 
jects populate the plane of ellipticity versus absolute B-magnitude 
(used to select the SAURON survey galaxies, see Paper II). 

A number of spectroscopic standard stars were observed to 
allow calibration of the stellar absorption line strengths to the 
Lick/IDS system (see Section l3!6l . Several photometric standards 
were also observed, giving an estimate of the absolute throughput 
and total spectroscopic response of the instrument and telescope 
over the three runs. 



Notes: (1) NGC number. (2) Hubble type fRC3: lde Vaucouleurs et alJl99lh . 
(3) Run number corresponding to Table fTI (4) Exposure time, in seconds. 
(5) Source for the seeing determination: either //Sr filter used or acquisition 
image (Acq.). (6) Seeing, full width at half maximum in arcseconds. 

2.3 Basic Data Reduction 

The data were reduced usin g the publicly available^ XOASIS 
software V6.2 <RousselllT993) developed at CRAL (Lyon). Pre- 
processing of the images included bias and dark subtraction. Dark 
frames were acquired for each exposure time used, including for 
calibrations and standard stars, to correct for non-linear 'hot' 
pixels on the CCD. The spectra were then extracted from the 
two-dimensional CCD image into a three dimensional data-cube by 
means of an instrument model, or 'mask'. This mask is essentially 
a map of the wavelengths falling on the different pixels for each 
lens, based on a numerical model of the instrument and telescope. 
The mask includes a model of the pupil for each lens, and simulates 
the curvature and chromatic variations of each spectrum on the 
CCD. Neighbouring spectra on the CCD are de-blended based 
on this model, and the spectra apertures are su mmed in th e 
cross-dispersion direction, applying optimal weights jHomell98d) . 

Although the extraction process imparts first-order wavelength 
information on the spectra from the instrument model, a more accu- 
rate wavelength calibration was performed in the conventional way 
using neon arc-lamp exposures: this removes higher-order varia- 
tions in the dispersion solution or inaccuracies in the mask. OASIS 
was mounted at the Cassegrain focus of the CFHT, and so flexures 

^ http://www-obs.univ-lyonl.fr/~oasis/download 



© RAS, MNRAS OQO.fTllsSI 



4 R.M. McDermid, et al. 



were accounted for by obtaining arc-lamp exposures bracketing the 
science exposures. Both neon exposures were used as a reference 
for the mask, giving the average position of the spectra on the CCD 
during the exposure. The dispersion solution itself is computed on 
a single arc-lamp exposure, generally using around 15 lines fitted 
with a second-order polynomial, with a mean dispersion of around 
a tenth of a pixel (~ 0.2 A). 

Flat-fielding of the data-cube was performed simultaneously 
in the three dimensions {x, y, A). Exposures of a continuum source 
were used for the spectral domain; exposures of the twilight sky 
were used to provide an illumination correction, or spatial flat-field. 
These two calibration exposures were combined to produce a single 
flat-field data-cube, which was divided through the science cubes. 

Cosmic ray hits affect neighbouring pixels on the CCD, which 
are effectively well-separated within the extracted data-cube, due 
to the staggered layout of spectra on the detector. This means that 
comparing spectra of neighbouring (on the sky) lenses (which 
should be rather similar due to PSF effects and the general nature of 
these objects) gives high contrast to cosmic ray events for a given 
spectral element. Pixels affected by cosmic rays were replaced by 
the median value of the neighbouring lenses at the corresponding 
wavelength. 

No sky subtraction was performed on these data since in all 
cases, the contribution from night sky lines was insignificant at the 
adopted wavelengths and within the small field of view centred on 
the bright galaxy nucleus. 

The spectral resolution can vary slightly within the field of 
view, which may introduce an irregular broadening of the lines 
when separate exposures are combined. We therefore homogenised 
the resolution to a common value. This was done by comparing 
twilight exposures with a high-spectral resolution solar reference 
spectrum'^ to measure the effective spectral resolution of each lens. 
These values were then smoothed with a circular spatial boxcar fil- 
ter of (X'S in radius to reduce noise. These 'maps' of spectral reso- 
lution were used to establish a representative instrumental broaden- 
ing, taken as 5.4 A (FWHM), the value for which 90% of the lenses 
(averaged over all available twilight frames of different nights and 
runs) were on or below that resolution. The twilight closest in time 
to the science exposure was used as a reference to smooth the data- 
cube with the appropriate Gaussian to obtain the adopted resolution. 

2.4 Fringe Pattern 

The data obtained during the three OASIS runs suffered from an 
interference pattern similar to fringing within the spectral range 
of our data. This added a periodic pattern with amplitude of about 
10 km s^^ to the mean absorption-line velocity fields derived from 
the spectra. The pattern was found to vary between exposures 
due to instrument flexure, and was therefore difficult to remove 
accurately. In principle, one would prefer to correct for this effect 
in the pre-extracted images, rather than on the extracted data-cubes, 
in which several adjacent pixels are combined at each wavelength. 
The interference pattern was not well understood during runs 1 
and 2, however, and so appropriate calibration frames (high-quality 
dome flats) were not obtained. For these runs, we were therefore 
required to make a correction on the extracted data-cubes. This 
was done by removing the average twilight spectrum (in which 
the effects of the fringing were largely randomised) from the 
twilight data-cube, leaving a residual fringe pattern, unique to 

^ http://bass2000.obspm.fr/solar_spect.php 
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Figure 2. Top left: mean velocity field of a twilight sky exposure derived 
using a solar reference spectrum and the pPXF method described in Sec- 
tion l3.1l A periodic pattern can be seen in the velocity field, resulting from 
a spatially-varying interference pattern in the spectra. Top right: mean ve- 
locity of the same twilight frame, corrected after extraction using a cor- 
rection frame constructed from a twilight exposure from a different night. 
This method was used to correct runs 1 and 2. Bottom left: mean velocity 
of the same twilight frame, corrected before extraction, using a median- 
combination of several dome-flat exposures. This method was used to cor- 
rect data from run 3. Bottom right: histogram of velocities from each field: 
black line = uncorrected field, blue line = dome-flat correction, magenta line 
= twilight correction. After the correction, the periodic pattern is removed, 
and the dispersion in the recovered velocities decreases from 7 km s~^ in 
the uncorrected case, to 4 km s~^ for both connection methods. (See text for 
details.) 



each lens. Dividing the science data-cubes by this frame removed 
much of the fringe pattern, although some signature was left due 
to flexures between the twilight and science exposures. For run 3, 
we obtained the appropriate calibration frames to allow the fringe 
correction to be made before extraction. A median combination of 
six dome flats was used to generate a correction frame. This frame 
was obtained by dividing each column of the master dome flat 
image by a high-order spline, fitted in such a way as to follow the 
complex light distribution, but not fit into the fringe pattern itself. 
Before extraction, we divided all data by the resulting correction 
frame to remove the fringe pattern. Small linear offsets between 
the correction image and the data image (due to flexure) were 
accounted for by comparing the associated arc-lamp exposures of 
the two frames, and shifting the correction frame accordingly. 

Fig. |2| demonstrates the effect of the fringing and fringe- 
correction on the mean stellar velocity field derived from a run 3 
twilight sky exposure (S/N > 100), which should show a flat ve- 
locity field. Velocity fields of the same twilight exposure are shown 
after correction with both methods. The standard deviation of ve- 
locity values decreases from 7 km s~^ in the uncorrected case to 
4 km s^^ with either correction. A residual gradient remains in the 
corrected velocity fields, with an amplitude of ~ 5 km s^^ across 
the field. This is within the fundamental uncertainties of our wave- 
length calibration, and below the typical measurement error of our 
galaxy observations (see Section lTTl . Although the correction us- 
ing dome flats is preferred, the solution using twilight exposures 
yields comparable results. 
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Table 3. Observed photometric standard stars. 



Name 


Run(#obs) 


Calibration 


HD093521 


2(3),3(4) 


STIS+Oke 


G191-B2B 


1(2),2(2) 


Model 


Feige 66 


1(1),3(1) 


FOS+Oke 


BD+33 2642 


1(2) 


FOS+Oke 


HZ 44 


1(1) 


STIS 



Notes: The second column indicates the ran number, and in parentheses, 
the number of observations during that ran. Only one star was observed 
per night. The third colu mn indicates the source of the referenc e calibration 
spectram, as described in iBohhn. Dickinson & Calzettil l200ll) . 



2.5 Flux Calibration 

During each night of the three observing runs, photometric standard 
stars were observed, primarily in order to calibrate chromatic vari- 
ations in throughput and instrument response, but also to give an 
absolute flux calibration for the data. The observed stars are listed 
in Table|3| and include several repeat measurements within and be- 
tween runs. A single spectrum was extracted from each reduced 
data-cube by integrating all spectra within a circular aperture cen- 
tred on the star To ensure that the flux star spectrum was representa- 
tive of the total flux in our wavelength range without simply adding 
noise, the PSF was measured from the reconstructed image of the 
star, and a summation of the flux within an aperture three times 
larger than the PSF FWHM was made, weighted by the Gaussian 
PSF model. 

The observed photometric standard stars were chosen to have 
as few spectral features as possible, thus making them suitable for 
calibrating residual spectral response variations. Over the relatively 
small spectral domain of OASIS, such variations are expected to 
be low-order However, imperfections in the calibration of our data 
(e.g. due to the fringe pattern) may give rise to fluctuations on scales 
of tens of Angstroms. Such residuals are difficult to detect, given 
the typically low spectral resolution of available spectrophotometric 
standard star reference tables (e.g. Oke, .1990). One solution is to 
compare with higher-resolution theoretical models. Such a model 
is available for the well-known star G191-B2B from Bohlin ( 2000), 
which provides flux measurements at varying wavelength intervals, 
from around 20 A at the red end of our wavelength domain (where 
the spectrum is varying only slowly) to < 1 A inside absorption 
features. 

Fig.|3|shows a comparison of this model with our four OAS I S 
observations before any flux calibration correction was applied. The 
model is interpolated to have the same spectral sampling as our 
observations, and is fitted to the observations using the penalized 
pixel-fitting method described in Section lTTI As well as applying a 
Gaussian velocity broadening (analogous to spectral resolution) and 
Doppler shift, the fitting process allows the inclusion of polynomial 
terms. The (average) fit shown in the top panel of Fig.|3|is achieved 
using a second-order multiplicative polynomial term. The bottom 
panel of Fig.|3|shows the residuals of the four individual fits to the 
OASIS observations, showing that the second-order polynomial is 
a good description of the relation between our observations and the 
theoretical model of G191-B2B. Deviations from the parabolic re- 
lation, although systematic across the two runs, are small, with a 
standard deviation of less than half a percent. This shows that struc- 
tures on intermediate wavelength scales which may be present in 
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Figure 3. Top panel: The four normalised OAS I S observations of the pho- 
tometric standard star G191-B2B (thin lines: magenta = run 1, green = run2) 
overplotted with the (normalised) average best-fit of the ,Bohlia (.2OO0i) theo- 
retical model (thick black line) using penalised pixel-fitting (pPXF, see Sec- 
tion l3.H . A multiplicative second-order polynomial term is included in the 
fit, as well as a varying shift and broadening of the model spectrum. Bottom 
panel: ratio of the OASIS observations and individual pPXF fits expressed 
as percentage deviations around the parabolic polynomial. The line colours 
indicate different runs as in the upper panel. Dotted horizontal lines indicate 
the (robust) standard deviation of the residuals with respect to zero, having 
an amplitude of 0.44%. 

our spectra due to e.g., filter throughput, are well removed by our 
spectral flat-fielding process. 

To evaluate whether there are global changes in the chromatic 
response of the observational setup, Fig.|4]compares the correction 
polynomials derived from all observed standard stars. In general, 
the curves agree well, with a maximum Icr spread of 2% at 
the spectra edges. Within a run, the scatter is significantly less 
(-^ 0.5%), however, and there appear to be some small systematic 
effects between runs. We therefore employ a single flux calibration 
curve per run, rather than a single correction for all data. 

We perform the flux calibration of our data compar- 
ing our observed stellar spectra w ith the calibrated spectra of 
iBohlin. Dickinson & Calzettil j200ll) . available from the "CAL- 
SPEC" database* of standard flux calibration spectra used for HST. 
Variations in the absolute flux calibration are as large as 20%. How- 
ever, to improve observing efficiency, no effort was made to obtain 
an accurate absolute photometric calibration by observing standard 
stars close (in time) to each galaxy pointing. 

2.6 Merging and Spatial Binning 

Two or more exposures were obtained for each galaxy (Table |2), 
each one offset by a small, non-integer number of lenslets to pro- 
vide oversampling and avoid systematic effects due to bad CCD 
regions. Multiple exposures were combined by first truncating to a 
common wavelength region and centring the spatial coordinates on 
the galaxy nucleus using images reconstructed from the data-cubes 
by summing flux in the spectral direction. Exposures were then re- 
normalised to account for transparency variations, and resampled 
from the original hexagonal spatial grid (set by the lens array) onto 
a common spatial grid of 0'.'2x0'.'2. Co-spatial spectra were com- 
bined via an optimal summing routine, taking into account the error 

^ http://www.stsci.edu/hst/observatory/cdbs/calspec.html 
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Figure 4. Normalised flux correction curves determined from each observed 
photometric standard star for all three observing runs (magenta = run 1, 
green = mn2, blue = run 3). From the stellar spectra, and from the compaii- 
son of G191-B2B shown in Fig.|3] a second-order polynomial was deemed 
sufficient to describe the chromatic response of the observational setup. 
Variations between the curves have a maximum Icr spread of 2% (with a 
range of 3%), although variations within a single ran are smaller (~ 0.5%). 
The systematic differences lead us to apply a single flux correction within a 
run, rather than a unique correction for all data. 

spectra which are propagated through the reduction. For the spec- 
troscopic standard stars, single spectra were extracted by summing 
all the data-cube spectra within a 2" circular aperture. 

At the stage of merging multiple pointings, a correction was 
made for differential refraction by the atmosphere, which causes 
the image to change position on the lens array as a function of 
wavelength. We applied the theoretical model of FilioDenko ( 1982), 
which uses atmospheric parameters to estimate the shift of the im- 
age as a function of wavelength and zenith angle. We adopted global 
values of 620 mb and 4°C for the air pressure and temperature, re- 
spectively, for all runs. Most observations were taken at airmass 
<: 1.3, ensuring that the typical relative shift due to atmospheric 
dispersion across the short wavelength range of our OAS I S obser- 
vations is < 0'.'2. No correction was made for atmospheric disper- 
sion for the standard stars, as the flux is combined over a relatively 
large aperture. 

It is necessary to analyse the galaxy spectra with a minimum 
signal-to-noise ratio (S/N) to ensure accurate and unbiased mea- 
surements. This was achieved using the binning method developed 
by CaDoellari & Cooin (2003), in which spectra are coadded start- 
ing from the highest S/N lenslet, and accreting additional lenslets 
closest to the current bin centroid. When the target S/N is reached, 
a new bin is begun. The resulting bin centroids are then used as 
the starting points for a centroidal Voronoi tessellation, generating 
compact, non-overlapping bins with little scatter around the im- 
posed minimum S/N. We conservatively bin all data to a minimum 
mean S/N of 60 per spectral resolution element for our analysis. 



2.7 Point Spread Function 

By comparing images constructed by integrating the flux in each 
spectrum with imaging from HST (where available), an estimate 
of the PSF can be made. We do this by parameterising the OAS I S 
PSF with two concentric circular Gaussians, and minimizing the 
difference between the HST image convolved with this PSF, and 
the reconstructed image. In practice, HST images in the appropriate 




flS 1-fl 1.5 2Si £5 

SAURON FWHM (areKtc) 



Figure 5. Comparison of OAS I S and SAURON PSF estimates for the galax- 
ies in Tablel2l The solid line indicates the 1 : 1 relation, showing that the PSFs 
of the OASIS observations are generally significantly smaller than those of 

SAURON. 

band are not available for all galaxies, in which case features such 
as dust may introduce small inaccuracies in our recovered values. 
For those galaxies with archival //5r images, the resulting PSF esti- 
mates, characterised by the FWHM of the best-fitting single Gaus- 
sian to our PSF model, are presented in column 6 of Table|2| Fig.|5| 
compares these values with the equivalent PSF measurements for 
the SAURON observations of these galaxies (Table 3 of Paper III). 

Although the spatial resolution of the SAURON and OASIS 
data sets are limited by atmospheric seeing conditions, the OAS I S 
data has the advantage of fully sampling even the best conditions 
on Mauna Kea within this wavelength range. For this reason, the 
OASIS data has generally better effective spatial resolution, and 
in some cases is significantly better, with more than a quarter of 
our sample having a PSF <. O/.'S. The median PSF of the OASIS 
observations (0'.'91) is almost half that of SAURON (l'.'7), with 
up to a factor three improvement in some cases. For four galax- 
ies the measured PSF of OASIS is similar to that of SAURON: 
namely NGC 1023, NGC 2699, NGC 4459 and NGC 4150. We nev- 
ertheless present the data in this paper, as the spectral domain of 
OAS I S includes several Lick iron indices which are not obtained 
with SAURON. 



3 DATA ANALYSIS 
3.1 Stellar Kinematics 

Stellar absorption-line kinematics were derived for each galaxy by 
directly fitting the spectra in pixel- space, using the penalised pixel - 
fitting method (hereafter pPXF) of ICappellari & EmseUeii] <2004 . 
as in Paper III. This method was chosen over Fourier-based methods 
due to its robustness to contamination by nebular emission lines 
(emission lines are simply excluded from the fitting process), which 
can be strong in our spectral domain around the central regions of 
early-type galaxies (e.g.. Paper V). 

Although robust to emission-line contamination, pixel-fitting 
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methods are sensitive to template mismatch e ffects, which can sig - 
nificantly bias the resulting kinematics (e.g.. iRix & Whitj|l992h . 
This problem was minimised by selecting an 'optimal template' 
at each iteration of the fitting process. For each trial set of kine- 
matic parameters, an optimal linear combination of appropriately 
convolved absorption spectra were fitted to the data. We use the 
same library of absorption spectra as in Paper III, which was taken 
from t he single-burst stellar population (SSP) models of Vazde^ 
jl999), with the addit ion of s everal individual stellar spectra with 
strong Mgb from the Ijone j il997h library (from which the SSP 
models are built) to compensate for the near-solar abundance ratio 
inherent in the SSP models. It was found that this library was able to 
reproduce the galaxy spectra very well, minimizing biases caused 
by template mismatch. 

After rebinning the wavelength axis of the data-cube to a log- 
arithmic scale (equivalent to pixels of constant velocity scale), we 
used the pPXF method to determine the line-of-sight velocity dis- 
tribution (herea fter LOSVD) parameterised by the truncated Gauss- 
Hermite series <van der Marel & Franxll993l:lGerharJl993h . Only 
the first four terms (mean velocity V; velocity dispersion a; and the 
higher-order terms /13 and /14, related to the skewness and kurtosis 
respectively) were measured, as these are generally sufficient to 
describe the significant features in early-type galaxy LOSVDs. Our 
data are generally of adequate quality to also measure fes and /le, 
which may be useful for some specific purposes (for example in the 
fitting of dynamical models), however they are not reported here. 

The penalty term of the pPXF method, quantified by the near- 
unity parameter, Ap, suppresses the large uncertainties inherent in 
measuring the higher-order terms /13 and /i4 when the LOSVD 
is undersampled by biasing the solution towards a simple Gaus- 
sian. In general, the objects in our sample have sufficiently large 
velocity dispersions such that this biasing is not important. Sev- 
eral objects, however, have velocity dispersions close to or be- 
lo w the instrumental dispersion. F ollowing the procedure described 
in ICappellari & EmsellemI (2004), we therefore optimised Ap via 
Monte Carlo simulations to suppress the measurement scatter to a 
level at which the true value (known in the simulation) lies well 
within the standard deviation of the biased value. This was found 
to be Ap — 0.5 for our adopted S/N of 60, differing slightly 
from the value used for the SAURON data in Paper III due to 
the different spectral resolution and wavelength coverage of the 
two instruments. From these same simulations, errors on the kine- 
matic parameters (V, a, /13, /14) for 5'/A'^ ~ 60 were found to be 
9 km s~^, 15 km s~^, 0.06, 0.06 respectively at a velocity disper- 
sion of UOkm s~^ and 12 km s"\ 12 km s"\ 0.04, 0.04 respec- 
tively at a velocity dispersion of 250 km s^^. 

Final errors on the derived galaxy kinematic parameters were 
determined via 100 Monte-Carlo realisations of each galaxy spec- 
trum, in which a representative spectrum of simulated white noise 
was added on top o/the noisy galaxy spectrum for each realisation. 
This approach, in contrast to adding noise to a noise-free broadened 
template spectrum, has the advantage that biases in the random er- 
ror due to template mismatch are included (although the systematic 
uncertainty caused by the mismatch is not). The difficulty in quan- 
tifying systematic uncertainties due to template mismatch is well 
known, and the reader is directed to Appendix B3 of Paper III for 
relevant discussion. 

Spectral regions which could possibly be affected by ionised 
gas emission were excluded from the fit by imposing 'windows' 
around the main features (H/3, [OIII] and [Nl]). The half-width of 
these windows was conservatively taken as 600 km s~^. In cases 
where the velocity dispersion of the detected gas agas was very high 
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Figure 6. Top: Comparison of OASIS and SAURON velocity dispersions 
measured on spectra integrated within equivalent circular apertures of 4" 
radius. Different coloured symbols indicate galaxies observed on different 
runs: magenta asterisk = Run 1, green diamonds = Run 2, blue triangles = 
Run 3. The dotted line indicates the unity relation, and the thick solid line is 
a robust straight-line fit to the data, with slope 0.97it0.02. Bottom: Differ- 
ence between the SAURON and OASIS velocity dispersion measurements, 
with quadratically summed errors. 



(taken as when iTgas > 400 km s~^), the windows were adjusted 
such that the half- width was given by 2 x (Xgas, and the fit was 
repeated. In practice, only NGC4486 (M87) exhibits such strong 
emission, and on the second iteration of pPXF, the stellar continuum 
fit was generally greatly improved. 



3.2 Comparison with SAURON stellar kinematics 

Due to the relatively small field of view of OASIS, the conven- 
tional comparison with literature values of velocity dispersion is not 
straight-forward, since these values are often integrated within, e.g., 
half an effective radius: generally much larger than the OAS I S field 
of view. Rather than introducing large and uncertain aperture cor- 
rections, the most direct comparison can be made with the SAURON 
measurements themselves, which were already shown to be consis- 
tent with independent studies in Paper III. 

To reduce the effects of atmospheric seeing on the dispersion 
values, we integrate the spectra from both instruments within a 
4" circular aperture, and derive the velocity dispersion from this 
high S/N spectrum using pPXF, fitting only for the mean velocity 
and velocity dispersion. Fig. |6| shows the resulting comparison of 
OASIS and SAURON velocity dispersion measurements {top panel) 
and the resulting scatter around the unity relation {bottom panel). 
The values are consistent between the two studies, with a mean off- 
set of 0.84±7.8 km s~^. The spectra used in this comparison are 
also used to compare the gas properties and line-strengths below. In 
Appendix Icl we present a comparison of the OASIS and SAURON 
spectra themselves. 



© RAS, MNRAS OQO.fTllsSI 



8 R.M. McDermid, et al. 



3.3 Kinemetry 

To quantify the structures found in our st ellar velocity fields, w e 
employ the so-called 'kinemetry' method of|Krainovic et al.V200f?). 
This techniqu e is an extension of the classic F ourier analysis of pho- 
tometry (e.g. Ijedrzeiewskil Il987t Ipranx. Il lineworth & HeckmanI 
[l989) . which can be used to describe the radial behaviour of two- 
dimensional kinematic maps. Annuli with increasing radius are 
considered, and for each annulus, a Fourier expansion is performed 
to fit the changing amplitude as a function of polar angle. 

By determining the flattened (elliptical) annulus along which 
the higher moments of the expansion are minimized, one can mea- 
sure the characteristic 'flattening' of the rotation field, equivalent 
to approximating the rotation to that of a thin disk. The amplitude 
of the first term, ki , gives a robust estimate of the peak rotation as 
a function of radius, equivalent to the circular velocity in the limit 
of a true thin disk. By also fitting the phase of the first moment 
of the expansion, a robust measurement of the kinematic position 
angle, PAkin, is obtained, defined here as the angular position of 
the peak rotation. This parameter traces the rotation axis, and is an 
important diagnostic of the kinematic structure of the object, since 
changes in the rot ation axis are a s ensitiv e indicator of deviations 
from axisvmmetrv lKrainovic et all i2005l) . 

Profiles of PAkin and ki are presented for all galaxies in 
Appendix [E] and are used to support out interpretation of the 
structures found in our velocity maps. For example, an abrupt 
change of PAkin combined with a minimum in ki indicates the 
presence of a kinematically decoupled component in the stellar 
velocity field that is counter-rotating to the outer parts, or has a 
different rotation axis. Detailed analysis of the full set of kinematic 
moments of the complete SAURON sample of E and SO galaxies, 
using both SAURON and OASIS data where available, will be 
given in a future paper of this series. 



3.4 Gas Properties 

The spectral region covered by our data includes a number of pos- 
sible nebular emission lines, namely H/3, [OIII], and [Nl]. It is nec- 
essary to separate these emission lines, where present, from the un- 
derlying stellar absorption continuum in order to obtain accurate 
line-strength indices, as well as to study the emission-line distri- 
bution and kinematics themselves. We do this using the method 
of Paper V, which includes a model of the emission-line spectrum 
as part of an optimal template fit. The stellar kinematics are fixed 
at the values determined by excluding the emission regions from 
pPXF as above. With the stellar kinematics fixed, the emission-line 
spectrum, modelled as a series of independent single Gaussians for 
each emission line, is included in the template library. The veloc- 
ity shift and broadening of the Gaussians are determined via non- 
linear optimisation, and the best-fitting non-negative linear combi- 
nation of stellar spectra and emission-line models is found at each 
iteration. This method has the advantage that both absorption-line 
templates and emission-line model are fitted simultaneously, avoid- 
ing the need for 'windowing' possible emission line regions. This 
makes optimal use of the information in the spectrum, resulting in 
accurate absorption and emission line measurements. Paper V gives 
a full discussion of this method. 

We adopt the approach of first fitting the [OIII] lines alone, and 
then repeating the fit including the H/3 and [Nl] lines in the emission 
spectrum, but with their kinematics fixed to the values measured 
from [OIII]. This approach was chosen since measuring the [OIII] 
lines is less sensitive to template mismatch, and furthermore, all the 



objects with detected H/3 emission show signs of [OIII] emission, 
which is usually stronger. For this reason, we use [OIII] as a robust 
indicator of emission, before considering that other lines may be 
present. In Paper V we discuss the merits and pitfalls of fitting H/3 
and [OIII] independently. Here we adopt the conservative approach 
of always assigning the kinematics derived from [OIII] to the H/3 
line. Although this may differ slightly from the true H/3 kinematics, 
in practice these differences are small compared to our measure- 
ment errors, and in any case do not significantly bias the fluxes we 
measure. 

The detection or non-detection of gas is based on the ratio 
of the amplitude of the fitted gas line and the scatter of the 
residual spectrum after subtraction of the best-fitting stellar and gas 
templates. By using the residuals of the template fit, this 'amplitude- 
to-noise' ratio (hereafter A/N) includes both random noise and 
systematic errors (e.g. due to template mismatch) in the noise esti- 
mate. The limiting values of A/N below which our measurements 
become strongly biased or uncertain were determined from simula- 
tions in Paper V. For the unconstrained fit of [OIII], emission lines 
with A/N ^ 4 are considered reliable. For H/3, the fitting process 
has less freedom, being constrained by the [OIII] kinematics. In 
this case, emission lines with A/N ^ 3 are considered reliable. 

The detection of [Nl] emission is complicated by the fact that 
the library of templates used to describe the stellar continuum are 
of near-solar abundance ratio. The effects of template mismatch are 
therefore particularly acute around the Mg 6 feature (A = 5177 A), 
slightly blueward of the [Nl] doublet. To incorporate this into our 
detection threshold, when computing the A/N for fitted [Nl] fea- 
tures, we base the residual noise estimate TV specifically on a wave- 
length region extending across both the [Nl] and Mg b features. Fol- 
lowing Paper V, we impose that A/N ^ 4 for this feature. In addi- 
tion, [Nl] is only considered if both H/3 and (by proxy) [OIII] emis- 
sion is also present. We further constrain the fit of [Nl] by imposing 
the same kinematics as H/3 and [OIII]. 

Based on these detection thresholds, we have shown in Paper V 
that it is possible to derive typical sensitivity limits, in terms of the 
equivalent width EW of the (Gaussian) emission line. Specifically, 
from eq. (1) of Paper V: 



EW = 



A/N X crobBV2^ 
S/N 



(1) 



where (Jobs is the effective dispersion of the fitted emission 
line including the instrumental broadening. For an emission line 
with intrinsic velocity broadening crint ~ 50 km s~^, a"obs ~ 
y^^~r^H = 144 km s-^= 2.4 A for [OIII] (2.3 A for H/3). 
For our imposed threshold of S/N — 60 for the continuum, this 
gives an estimate for the limiting EW of [OIII] of 0.40 A for the 
chosen limit of A/N ^ 4 (0.29 A for H/3). For given A/N thresh- 
olds, the limiting EW for OAS I S compared with SAURON (within 
the same aperture) scales as the ratio of the instrumental broadening 
and S/N, such that 



EWa^,, ^ Ooas ^ S/Ns, 
SWsau (Tsau S / No. 



1.25 X 



S/N,, 
S/No. 



(2) 



In the central few arcseconds of our galaxies, the SAURON data 
are generally not spatially binned, and have typical S/N far in ex- 
cess of the binning threshold. The OAS I S lenses are smaller on the 
sky than those of SAURON, with more than a factor 10 less col- 
lecting area. For this reason, the OASIS spectra are usually close 
to or below the binning threshold in the same regions. The sensi- 
tivity in these regions is consequently less with OAS I S than with 
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SAURON, although the effective spatial resolution is higher (see also 
AppendixlDl. 

Errors on the emission line parameters are estimated via 100 
Monte Carlo realisations of the full analysis process. The galaxy 
spectra (including the emission) with added noise used in the sim- 
ulations to compute the kinematic errors (see Section im are used 
as input to the gas-fitting method. In this way, we propagate errors 
in the template-fitting process, as well as the emission line fitting 
itself. The amplitude of the errors varies with the A/N of the emis- 
sion line, which can vary significantly within and between galaxies. 
For A/N — 60, typical errors in the derived flux, velocity and dis- 
persion are 3%, 2% and 10% respectively. For A/N = 10, these 
values increase to 10%, 4% and 25% respectively. 

3.5 Comparison with SAURON gas measurements 

Given the complexities of separating often faint emission lines from 
the stellar continuum, and the high-quality data required to make re- 
liable measurements, making a meaningful comparison of our mea- 
sured gas properties with literature values is a difficult task. In Pa- 
per V we established that there is a satisfactory agreement between 
SAURO N gas measurements and those of the Palomar spectroscopic 
survey ilia Filippenko & Sargent.l995>) . given the differences in 
techniques and data quality. We therefore again adopt the SAURON 
measurements as our reference, since the data quality are similar to 
our observations, and the apertures used can be accurately matched. 

Fig.Qpresents a comparison of the [OlII]/H/3 ratio measured 
with SAURON and OASIS. This demonstrates our ability to mea- 
sure emission line fluxes, independent of the absolute flux calibra- 
tion of the two data sets, which are both approximate. This compar- 
ison shows reasonable agreement, with a mean SAURON:OASIS 
ratio of 1.19 ± 0.44. There is some evidence of a systematic trend, 
with a fitted slope of 1.08 ± 0.03 from the error-weighted line fit. 
Fig.|8|compares the velocity dispersion of the [OIII] doublet. Again, 
there is reasonable agreement, with a mean ratio of 0.91±0.24, and 
a fitted slope of 0.99 ± 0.05 from the error-weighted line fit. 



3.6 Line Strengths 

The spectral range of our OASIS data contains a number of 
absorption features which can be used as diagnostic tools to 
determine the distribution of stellar populations within a galaxy, 
in terms of the luminosity- weighted population parameters: age, 
metallicity and element abundance ratio. We quantify the depth of 
these absorption features through the use o f atomic line-strength 
indices, taken from the Lick/IDS sys tem jBursteinetailligS^ 
IWorthevetai]|l994 iTraeer et al1ll998h . the bandpass definitions 
of which are given in Table |4| 

Three absorption features in our spectral range, namely H/3, 
Fe50I5 and Mgfe, can be significantly altered due to 'infilling' by 
possible H/3, [OIII] and [Nl] emission features, respectively. Before 
measuring the absorption-line strengths, we therefore subtract the 
emission spectrum as derived in Section l7!4l from the original data. 
Finally, the absorptio n-line strengths are calibrated to the classical 
Lick/IDS system (e.g. lTrager et alJl998^ . Errors are estimated via 
100 Monte Carlo simulations of our complete line-strength analy- 
sis, using the output 'cleaned' spectrum from the Monte Carlo sim- 
ulations of the emission line analysis (see Sections |3. II and Ua\ . In 
this way, we include the uncertainty in our emission line correction 
and stellar kinematics directly. For S/N ~ 60, typical errors for the 
indices H/3, Fe5015, Mg b, Fe5270, Fe5335 and Fe5406 are 0.15 A, 
0.3 A, 0.15 A, 0.15 A, 0.17 A, and 0.13 A respectively. 




Figure 7. Top: Comparison of OASIS and SAURON [OIIIJ/H/3 measure- 
ments within a circular aperture of 4" radius. The dotted line shows the 
unity relation, and the thick black line shows a straight-line fit weighted by 
the errors in both axes. Bottom: Ratio of the SAURON and OASIS mea- 
surements, with quadratically summed errors. Different coloured symbols 
indicate galaxies observed on different runs as in Fig.m 
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Figure 8. Top: Comparison of OASIS and SAURON [OIII] velocity disper- 
sion measurements within a circular aperture of 4" radius. The overplotted 
lines ai'e as in Fig.Q Bottom: Ratio of the SAURON and OASIS measure- 
ments, with quadratically combined errors. Different coloured symbols in- 
dicate galaxies observed on different runs as in Fig.l6l 
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Table 4. Bandpass definitions of Lick indices within our OAS I S spectral 
range. 



Index 


Units 


Central bandpass 


Pseudo-continuum 


H/3 


A 


4847.875 - 4876.625 


4827.875 


- 4847.875 








4876.625 


-4891.625 


Fe5015 


A 


4977.750 - 5054.000 


4946.500 


- 4977.750 








5054.000 


- 5065.250 


Mgb 


A 


5160.125 -5192.625 


5142.625 


-5161.375 








5191.375 


- 5206.375 


Fe5270 


A 


5245.650-5285.650 


5233.150 


-5248.150 








5285.650 


-5318.150 


Fe5335 


A 


5312.125 -5352.125 


5304.625 


- 5315.875 








5353.375 


- 5363.375 


Fe5406 


A 


5387.500-5415.000 


5376.250 


- 5387.500 








5415.000 


- 5425.000 



Below we describe the steps taken to calibrate our line-strength 
measurements to the Lick system. To account for the difference in 
spectral resolution of our observations and the Lick/IDS system, 
and to facilitate comparison with other authors, we convolve our 
spectra using a wavelength-dependent Gaussian kernel based on a 
linear interpolation of the resolution values given in the appendix 
of .Worthev & OttavianL L1997.) . 

3. 6. 1 Velocity Broadening Correction 

The absorption lines of integrated galaxy spectra are broadened 
by the collective velocity distribution of the stars along the line 
of sight, usually approximated by the velocity dispersion, a. The 
effect of this broadening is generally to weaken an index by lower- 
ing the surrounding continuum bands, as well as contaminating the 
central pass-band with contributions from neighbouring absorption 
features. The conventional correction for this effect involves 
broadening a collection of observed stellar spectra and measuring 
the Lick indices as a function of a only. The size of this correction 
can vary strongly as a function of a and, for some indices, with the 
intrinsic strength of the absorption feature. In addition, deviations 
from a Gaussian velocity distribution can also introduce a small 
systematic error following this approach tKuntschner,.,2004) . 

We therefore refine this method by using the optimal template 
derived for the stellar kinematics to determine the velocity broad- 
ening correction of each spectrum (as in Paper VI). The correction 
is given by comparing the index values of the unbroadened optimal 
template with those of the template broadened by the derived 
velocity distribution, including higher-order moments. The optimal 
templates derived by pPXF provide very good representations of 
the galaxy spectra in our sample, thus the method accounts directly 
for dependence on intrinsic line strength, as well as the detailed 
shape of the velocity distribution, without relying on average 
corrections determined from a range of stellar or model templates. 

3.6.2 Lick Offsets 

The Lick system is based on non-flux-calibrated spectra, and it is 
therefore necessary to account for offsets in the measured line in- 
dices caused by small differences in the local continuum shape of 
our spectra when compared with those of the (non-flux calibrated) 
Lick system. This is done by observing a number of stars in com- 



Table 5. Estimate of Lick offsets derived from all stars. 



Index 


Offset 


H/3 


-0.10 ±0.03 A 


Fe5015 


0.45 ± 0.07 A 


Mgfe 


-0.02 ± 0.03 A 


Fe5270 


-0.04 ± 0.03 A 


Fe5335 


-0.03 ± 0.04 A 


Fe5406 


0.05 ± 0.04 A 



Note: For each index, the biweight mean and dispersion estimate was used to 
reduce the influence of outlying values. The error on the mean offset derived 
from the observations is given as the dispersion scaled by 1/ \/~N. 

mon with the Lick librar y, and comparing the m easured indices with 
those in the appendix of lWorthev et alJil994) . 

During the OASIS observing campaign, a total of 31 Lick 
standard stars were observed, including 8 repeat observations be- 
tween one run or more, and 3 repeat observations within runs. In 
total, there were 45 observations of stars in common with the Lick 
library. Fig.|9|presents the differences between the reference Lick 
values and our observations, using all measured stars (repeat mea- 
surements are included as separate points). Unfortunately there are 
too few repeat observations to fully evaluate night-to-night and run- 
to-run variations. However, the maximum difference between ob- 
servations of the same star are consistent with the scatter of all 
measurements. From the distribution of points in Fig.|9| the offsets 
measured from each run individually are self-consistent, suggest- 
ing minimal run-to-run biases. The final offsets are given in Table 
|5| which have been added to the measured indices presented here 
unless otherwise stated. 

In general, the determined offsets are small, and consistent 
within the Icr errors of those found in Paper VI, where many more 
stars were used. The exception to this is Fe5015, which shows an 
offset larger than expected from SAURON, consistent only at the 
2a level. For all indices, a number of stars are clear outliers of 
the distribution. Often these are M-type stars, where small errors 
in the velocity measurement (required to place the bandpasses at 
the appropriate rest wavelengths) can give rise to large changes in 
the measured index. The mean offset is determined with a biweight 
estimator ( Hoaalin. Mo steller & Tukev. .l983.') . and so is robust to 
these deviant points. 



3.7 Comparison with SAURON line strengths 

As for the comparisons with SAURON data presented in Sections 
13.21 and 13.51 the line-strengths are best compared when derived 
from spectra within the same circular aperture from both data sets. 
Fig. llOl presents a comparison of OAS I S and SAURON line-strength 
measurements within the same 4" aperture used before. Due to 
truncation of the SAURON wavelength range by the instrument 
filter, only three lines can be measured within this aperture for all 
galaxies: H/9, Mgb and Fe5015. The measurements are made after 
the emission lines have been removed as described in Section l7!4l 
Given the sensitivity of line-strength measurements to systematic 
errors, as well as the inherent difficulty of separating the emission 
lines from the stellar absorption features, the agreement between 
the two data sets is remarkably tight. For H/3 and Mg b, the 
(biweight) mean difference and standard error on the mean are 
-0.09 ± 0.03 A and 0.05 ± 0.04 A respectively. For Fe5015, the 
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Figure 9. Measured differences between Lick index reference values iWortliev et alll994 and OAS I S observations for the six indices in our spectral range. 
Different coloured symbols indicate stars observed on different runs as in Fig.|^ Repeat measurements of the same star are included as separate points. The 
biweight mean and robust dispersion estimates based on all measurements are indicated by the dashed and dotted lines respectively. Values of the biweight 
mean and the error on this mean are given in Tablelsl 




2 3 



5-^ ' '- 

Q ^ 

J - ' 
fir. 




X! - 

5 r 

^L.„j 





4.0 



4.5 5.0 
Ee50:5L 



5,5 
(A) 



6.0 6.5 




2 3 



Figure 10. Comparison of SAURON and OASIS line-strength measurements. Different coloured symbols indicate galaxies observed on different runs as 
in Fig.|6] Dotted lines show the identity relation; solid lines indicate robust straight-line fits, from which slopes are derived. Since both data sets are flux 
calibrated, no Lick offsets are applied to the values in this comparison. For both H/3 and Mg b, the mean differences are consistent with zero, with the fitted 
slopes consistent with no measurable systematic trends (1.0 it 0.04 and 0.95 ± 0.05 respectively). The agreement is less good for Fe5015, which gives a fitted 
slope of 1.0 ± 0.09, but with a constant mean offset of -0.36it0.05A. 
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agreement is less good, with a mean difference of —0.36 ± 0.05 A, 
and evidence for some systematic behaviour. 

3.8 Age, Metallicity and Abundance Ratio 

The purpose of measuring the line-strength indices is ultimately to 
understand the distribution of stellar ages, metallicities and element 
abundances in these galaxies. Although certain features can be more 
sensitive to these quantities than others, there is always a certain 
degree of degeneracy between the population parameters and sin- 
gle absorption features. By using a combination of line-strength in- 
dices, it is possible to constrain stellar population models, which in 
turn allow us to interpret the distribution of line-strengths as one of 
age, metallicity and element abundance. 

In order to determine these physical quantities, 
we use t he single-burst stella r population models of 
iThomas. Mar aston & BendeJ i2003h . We interpolate the origi- 
nal grid of models to produce a cube of 40 x 40 x 40 = 64000 
individual models, spanning [Z/H] = —0.33 to +0.67 (in steps of 
0.026 dex), age = 0.1 to 15 Gyr (logarithmically, in steps of 0.056 
dex), and [o/Fe] = -0.2 to 0.5 (in steps of 0.018 dex). We locate 
the model which lies closest to our six measured Lick indices 
simultaneously for each bin in our data using the technique, sim- 
ilar to the method described in Proctor. Forbes. & Beaslev (2004). 
We estimate the confidence levels on the derived parameters using 
the Ax^ values, which were checked to be consistent with those 
found from our Monte-Carlo simulations. For measurements which 
lie outside of the model population parameter-space, values on the 
closest boundary are chosen. 

Fig. ll ll shows an example of the Ax^ contours obtained with 
this method, indicating the minimum and confidence levels. This 
shows that the typical Icr uncertainties we obtain on the parameters 
of metallicity, abundance ratio and age are of the order 0. 1 dex. The 
effect of the age-metallicity degeneracy (e.g., Worthev et al. 1994) 
can clearly be seen as a tilt in the contours of the right-most panel 
ofFig.EI] 

We note here that the approach of using SSP models is intrin- 
sically limited to providing only luminosity-weighted properties. 
Moreover, this weighting does not affect all indices equally: for 
example, H/3 is strongly enhanced at younger ages (0.5 - 5 Gyr), 
while other diagnostic features like Mg and Fe are less so. As a re- 
sult, it is not certain that the different SSP parameters derived (i.e. 
age, [Z/H], [of/Fe]) are weighted in the same way from the different 
populations that may be present. One must therefore be cautious in 
interpreting the maps of population parameters, since the different 
maps may be sensitive to different populations. 



4 RESULTS 

Maps for all measured quantities (stellar kinematics, ionised gas 
fluxes and kinematics. Lick indices and stellar population param- 
eters) are presented in Fig. ^| To help give an impression of our 
measurement errors, which can be difficult to infer only from the 
scatter in the parameter maps themselves, we present a simulated 
major-axis 'slit' of NGC 3489 in Appendix |^ showing each of 
the measured and derived quantities and their individual error bars. 
There is a wealth of structure in the parameter maps, much of which 
is unseen or poorly resolved in our SAURON data. For a compre- 
hensive description of the global properties and bibliographic refer- 
ences of the galaxies presented here, the reader is directed to Papers 
III, V and VI of this series. In this section, we present an overview of 



the central structures found in the various parameter maps, consid- 
ering the stellar kinematics, gas properties, and stellar populations 
in turn. 



4.1 Stellar Kinematics 

The variety of stellar motions observed in the SAURON sample 
can be broadly separated into two categories: objects that exhibit 
rotation (and possible counter-rotation) around a single axis; and 
objects that have a varying axis of rotation. Based on the large (kpc) 
scales probed by the SAURON observations, Emsellem et al. (Paper 
IX, in preparation) refine this picture, and quantify the behaviour in 
terms of the specific angular momentum of the galaxy, finding that 
fast-rotating galaxies exhibit a rotation axis which generally does 
not vary; whereas slowly-rotating objects tend to show misaligned 
components. The differentiation between fast and slow rotating 
galaxies is visually apparent from the velocity maps in Paper 
III. On the scales probed by OASIS, we find that a number of 
objects, while kinematically well-aligned on large-scales, may still 
exhibit misaligned central components. For this reason, we divide 
the following qualitative discussion between galaxies exhibiting 
misaligned and aligned rotation within the OASIS maps, rather 
than between fast and slow rotating objects. 

In addition to the maps, we present results from our kinemet- 
ric analysis in Appendix |n| giving for each galaxy the kinematic 
position angle and maximum absolute rotation as a function of 
galactocentric radius. We use the kinemetric profiles of kinematic 
position angle (PAkin) to help determine the velocity fields which 
show misaligned and aligned rotation. The kinemetric profiles are a 
useful tool to help quantify structures which can be approximately 
inferred from the maps, although, as with photometric profiles, one 
must also bear in mind the effects of seeing and spatial sampling 
when interpreting them. When the velocity gradient is low, the error 
in PAkin is also larger. With these caveats in mind, we consider 
galaxies as being aligned if the variations in PAkin are around 10° 
or less within the OASIS field of view. For APAkin> 10°, the 
velocity field is considered misaligned, with strong misalignments 
(APAkin J> 30°) indicating a decoupled component. A more 
complete definition of kinematic structures found in our SAURON 
and OASIS maps will be given in a future paper of this series. 



4.1.1 Misaligned Rotation 

Several of the velocity maps presented in Fig. 1121 show complex 
velocity fields, where the axis of rotation changes as a function 
of radius. This is illustrated by the variation in the PAkin profiles 
in Appendix IHl which range from smooth variations of a few de- 
grees, to an abrupt change of many tens of degrees. In addition, 
some galaxies rotate around an axis which shows little variation in 
orientation, but which is misaligned with respect to the local pho- 
tometric axes. Such features may indicate the presence of a distinct 
dynamical sub-component, such as a bar, or may be indicative of a 
triaxial potential (e.g. .Statler.l991.) . 

Differentiating between these alternatives is not always 
straightforward, especially when probing the limits of spatial and 
spectral resolution. The orientation of the system also plays an 
important role in determining the projected quantities we observe. 
With this in mind, we identify several manifestations of kinematic 
misalignment present in our data, varying from subtle 'twists' in 
the zero-velocity contour (ZVC), to spatially distinct components 
with an apparently different orbital composition from that of 
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Figure 11. Example of Ax^ space for finding the best-fit model to the lines strengths of NGC 4564, as measured from the 4" aperture spectrum described in 
previous sections. The best-fit model is found by minimizing the difference between the model predictions and the six Lick indices from our data, measured 
using a approach. The three panels show different projections of the three-dimensional solution space (age, metallicity [Z/H], abundance ratio [ct/Fe]), with 
one of the parameters fixed at the best-fit value. The thick contours indicate the \a, 2cr and 3(t confidence levels, based on Ax^ for two degrees of freedom. 
We note that in the plane of age and metallicity (far right panel), the effect of the age-metallicity degeneracy can be seen as the tilt of the contours with respect 
to the plot axes. 



the rest of the galaxy. Moreover, some objects appear regularly 
rotating in the OAS I S field of view, but are misaligned at the larger 
scales observed with SAURON. These objects are included in this 
discussion, as the OAS I S observations can provide a more detailed 
perspective on the nature of the central component. 

The following galaxies show mild twists (10° < APAkin 
< 20°) in their ZVC within the OASIS field of view: NGC 1023, 
NGC 3489, NGC 4459 and NGC 4552. In the case of NGC 1023, 
the axis of rotation appears to change gradually as a function of ra- 
dius even to the edge of the SAURON field (Paper III). In the other 
cases, the twist is more localised to the central regions. The strongly 
barred galaxies NGC 2699 and NGC 4262 show a single axis of ro- 
tation that is constantly misaligned with the local photometric prin- 
ciple axes, though the central kinematic axes align quite well with 
the large scale photometry. 

More pronounced twisting of the ZVC is found in both 
NGC 2768 (APAkin~ 40°) and NGC 5982 (APAkin~ 30°). The 
central regions of these objects show a change in rotation axis 
from minor- to major-axis. This change is clearer in NGC 5982, 
where the central component is quite distinct from the outer 
body, showing a peak in the profile of maximum rotation (ki) in 
Appendix |n| In NGC 2768, the rotation increases smoothly as 
a function of radius, with no distinct rotation component in the 
centre. The velocity dispersion, however, shows a clear dip within 
the region of the twisted ZVC, showing that the central kinematics 
are different from the outer parts. 

These two objects mark something of a transition between 
twisted velocity fields and decoupled components. In the case of 
mild twists, the point at which the galaxy exhibits a dynamically 
distinct component is almost impossible to identify with any degree 
of certainty, at least from our data alone. In NGC 2768, the twist is 
stronger and localised within 1" of the centre. The link with a sharp 
dip in the velocity dispersion gives the strong impression that the 
orbital distribution is indeed different from the rest of the galaxy. 
In NGC 5982, the case for a distinct central component is stronger, 
given that there is a peak in the rotation curve, coincident with 
some lowering of the velocity dispersion, suggesting a possible 
dynam ically cold component. For this specific galaxy, Statler 
(l991) highlighted the important issue that a triaxial system with a 
smooth distribution function, when seen in projection, can naturally 



give rise to a velocity field very similar to NGC 5982 (see also 
lOosterloo. Balcells & Carte3 'l994'). Whether or not components 
are actually decoupled in phase-space r equir es the application of 
dynam ical models (such as Icappellari eiraDl200il : IVerolme et alj 
i2003t:iStatIer et all2004h . which is beyond the scope of this paper. 

Galaxies which have clearly distinguished misaligned compo- 
nents (i.e. APAkinS^ 40°), which we will term kinematically decou- 
pled components (hereafter KDCs), include NGC 3032, NGC 4382, 
NGC 4621, and NGC 5198. Of these, the misaligned components in 
NGC 3032 and NGC 4382 were not resolved in our SAURON data 
(Paper III) and are reported here for the first time. Photometrically, 
NGC 4382 exhibits a peculiar core region, showin g a central d rop in 
surface brightness that is not explained by dust i Lauer et all2005) . 
Figure fT2l shows this galaxy has a residual rectangular central struc- 
ture in the unsharp-masked /fST image, which is slightly misaligned 
with the outer isophotes (by ~ 10°). This boxy structure is spatially 
coincident with the KDC and oriented in a similar direction, sug- 
gesting that the two components are directly related. 

NGC 4621 is also visible a s a KDC in our data. It was first dis- 
covered bv lWemli et alj i2002h using calcium II triplet data taken 
with OASIS (CFHT) using adaptive optics. NGC 4150 shows evi- 
dence of a KDC in the SAURON data of Paper III, but the poor see- 
ing conditions during our OAS I S observations preclude an accurate 
measurement of this component in our data. With the exception of 
NGC 5198, each of these apparently small KDCs reside in galaxies 
with otherwise regular, aligned rotation at larger radii. The flattest 
of these objects, NGC 4621, shows clear evidence that the KDC is 
embedded in a dynamically cold disk, visible in the residuals of the 
unsha rp-masked HST image in Fig. ll2l (see also iKra i novic & Jaff j 
l2004h . 

NGC 3414, NGC 3608, NGC 5813, and NGC 5831 contain 
large KDCs, clearly visible in the SAURON observations (Paper 
III). These decoupled components fill the OAS I S field of view, and 
generally show regular rotation, with some possible twisting of the 
ZVC as the field encroaches on the outer component beyond the 
OAS I S field of view. Both NGC 3414 and NGC 5813 exhibit strong 
rotation and anti-correlated /13 fields, as well as depressed velocity 
dispersion where the rotation is highest, indicating that these KDCs 
are disk-like. 

The variety of objects with misaligned rotation emphasises 
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the fact that such behaviour can result from a number of different 
intrinsic orbital structures and possible sub-components. In some 
cases, a large-scale bar structure gives rise to the twisted velocity 
field; in others, the decoupled component takes the form of a disk 
embedded in a slowly-rotating spheroid; and in some galaxies, the 
decoupled component cannot be resolved even at sub-arcsecond 
scales, making the intrinsic structure unclear. Grouping these 
objects together belies the different origins of the misalignment 
(see also|5}, and also highlights the difficulty in separating galaxies 
containing distinct dynamically decoupled components from those 
having a complex, but smoothly varying orbital structure. With 
the application of new dynamical modelling techniques to study 
galaxies with twists and KDCs, we hope to obtain a clearer under- 
stand ing of these objects' orbital structure (e.g. Ivan de Ven et alJ 
|2004. van den Bosch et al., in preparation). 



4.1.2 Aligned Rotation 

The remaining objects in our sample show little or no evidence of 
misaligned stellar rotation in the OASIS field of view. The two 
giant elliptical galaxies, NGC4486 and NGC5846, are included 
here, as they show almost no rotation. The strongly barred galaxy, 
NGC 3384, is also included here, as the rotation is very symmetric 
around the local photometric minor axis. Due to the morphology of 
the bar, the photometric PA shows some rapid changes as a func- 
tion of radius, but becomes aligned with the central PA again in the 
outer parts ( Busarello et al. 1996). Within our OAS I S and SAURON 
fields, the rotation axis remains quite constant. In the OAS I S data, 
this galaxy exhibits evidence of a stellar disk, with strong rotation 
and anti-correlated hj, fields, and a remarkable depressed velocity 
dispersion along the major axis. This is consistent with the disk-like 
residuals in the unsharp-masked HST image. 

The remaining objects (NGC 2549, NGC 2695, NGC 2974, 
NGC 3379, NGC 4473, NGC 4526, NGC 4564 and NGC 5845) ap- 
pear consistent with rotation around a unique axis ( APAkin 10°), 
even on the larger scales probed by SAURON. The galaxies which 
show the most flattened isophotes within the OAS I S field - namely 
NGC 4473, NGC 4526 and NGC 5845 - show strong evidence of 
harbouring central stellar disks. This evidence comes from disk-like 
residuals in the unsharp-masked HST image, strong rotation fields 
with 'pinched' isovelocity contours, significant h-j, amplitude anti- 
correlated with the velocity, and lowered velocity dispersion in the 
region occupied by the disk. In this last characteristic, NGC 4473 
is peculiar, showing a rising velocity dispersion along the major- 
axis. Detailed dynamical modelling of this galaxy indicates that it 
is composed of two flattened, almost equal-mass, components ro- 
tating in opposite directions iCaDoellari & McDermid 2005). The 
central drop in the dispersion indicates where the contribution from 
the counter-rotating component becomes less significant. 



4.2 Ionised Gas 

In Paper V we have shown that many of the galaxies in our 
sample contain large quantities of ionised gas, spanning a range 
of distributions from rotating disks to irregular filaments. Here 
we discuss the incidence of ionised gas within the central regions 
of our sample galaxies observed with OASIS, and describe the 
various distributions and kinematics we find. 



Table 6. Summary of central dust and ionised gas properties observed with 

OASIS. 



NGC 


[OIII] 


H/3 


[NI] 


Dust 


(1) 


(2) 


(3) 


(4) 


(5) 


1023 


trace 


no 


no 


trace (filaments) 


2549 


yes 


trace 


no 


no 


2695 


no 


no 


no 


no image 


2699 


no 


no 


no 


no 


2768 


yes 


yes 


yes 


spiral 


2974 


yes 


yes 


yes 


spiral 


3032 


yes 


yes 


yes 


spiral 


3379 


yes 


yes 


no 


ring 


3384 


no 


no 


no 


no 


3414 


yes 


yes 


yes 


trace (filaments/spiral) 


3489 


yes 


yes 


no 


spiral 


3608 


yes 


trace 


no 


no 


4150 


yes 


yes 


no 


spiral 


4262 


yes 


yes 


no 


trace (spiral) 


4382 


no 


no 


no 


no 


4459 


ti'ace 


trace 


no 


spiral 


4473 


no 


no 


no 


no 


4486 


yes 


yes 


yes 


trace (filaments) 


4526 


yes 


yes 


no 


edge-on dust lanes/spiral 


4552 


ves 


trace 


trace 


filaments 


4564 


no 


no 


no 


no 


4621 


trace 


no 


no 


no 


5198 


yes 


yes 


no 


trace (filaments) 


5813 


yes 


yes 


no 


filaments 


5831 


no 


no 


no 


no 


5845 


no 


no 


no 


central disk 


5846 


yes 


yes 


yes 


filaments 


5982 


yes 


trace 


no 


no 



Notes: (1) NGC number. (2) Presence of [OIII] emission. (3) Presence of H/3 
emission. (4) Presence of [NI] emission. (5) Presence of visible dust features 
in the unsharp-masked HST image, with some qualitative description of the 
dust morphology. The F555W filter was used wherever available. For all 
columns, the term 'trace' implies that features were found close to the limit 
of detection. 



4.2.1 Incidence of Ionised Gas 

The incidence of ionised gas in our sample is summarised in Table 
|6| We find that 14 of the 28 galaxies in our sample have clear evi- 
dence of both H/3 and [OIII] emission in their central regions, and a 
further four objects with only [OIII] clearly detected. Three galaxies 
exhibit faint [OIII] close to the detection limit (marked as 'trace' in 
Table|S), with either undetected or similarly faint H/3. Seven galax- 
ies show evidence of [NI] emission, which is generally confined to 
the central regions. 

This incidence of gas is comparable to that found with 
SAURON, although there are a number of galaxies for which 
SAURON shows stronger detection of emission, or shows clearer 
structures where the OASIS data appear patchy or even show no 
emission at all. This is simply a reflection of the sensitivity lim- 
its of the two data sets, arising mainly from the generally lower 
S/N of the OASIS spectra compared to the SAURON spectra, as 
discussed in Section l7!4l Appendix|D|quantifies this further, show- 
ing for the apparently discrepant cases the different 5'/A'^ values of 
the SAURON and OAS I S data sets as a function of circular radius, 
and how this translates into equivalent width limits for the two in- 
struments. The [OIII] emission detected by SAURON falls below the 
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NGC1023 




Figure 12. The following figures show OAS I S maps of (from top to bottom, left to right): ^>.vf column - reconstructed intensity, stellar velocity, stellar velocity 
dispersion, Gauss-Hermite coefficients and /14, second column - [OIII] emission-line equivalent width, [OIII] gas velocity, [OIII] gas velocity dispersion, 
H/3 emission-line equivalent width, [OIII]/H/3 emission line ratio, third column - unsharp-masked HST image (when available), Fe5015, Fe5270, Fe5335, 
Fe5406 Lick absorption-line indices, fourth column - H/3, and Mg b Lick absorption-line indices, luminosity-weighted stellar age, metallicity and abundance 
ratio. An arrow beside the figure title indicates north-east. The HST images are given as the ratio of the oiiginal image and the image convolved by a Gaussian 
kernel of 10 pixels. The overplotted isophotes are in half-magnitude steps. Grey bins in the emission-hne maps indicate non-detections. 
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Figure 12. continued 
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Figure 12. continued 
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Figure 12. continued 
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Figure 12. continued 
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Figure 12. continued 
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Figure 12. continued 
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Figure 12. continued 
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Figure 12. continued 
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Figure 12. continued 
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NGC3608 Sl- 




Figure 12. continued 
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Figure 12. continued 
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Figure 12. continued 
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Figure 12. continued 
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Figure 12. continued 
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Figure 12. continued 
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NGC44S6 Sl- 




Figure 12. continued 
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Figure 12. continued 
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Figure 12. continued 
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Figure 12. continued 
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Figure 12. continued 
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Figure 12. continued 
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Figure 12. continued 
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Figure 12. continued 
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Figure 12. continued 
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Figure 13. Central spectrum of NGC 4486, showing double-peaked emis- 
sion lines. This is indicative of unresolved gas rotating rapidly (~ 
±500 km s~^) in a nuclear disk sun'ounding the central black hole. 

equivalent-width limit of our OAS I S data, thus accounting for the 
difference in detection. 



4.2.2 Ionised Gas Kinematics: Rotation 

The mean motion of the ionised gas is characterised by three 
main morphologies: regular circular rotation, consistent with a disk; 
warped or twisted rotation; and irregular motion. We discuss each 
of these in turn. 

The following galaxies show gas kinematics consistent with 
disks: NGC 3379, NGC 3608, NGC 4486 (centre), NGC 4526, 
NGC4552, NGC5198, and NGC5813 (centre). Of these, 
NGC 4486 is worth a special mention. Firstly, the PSF of these 
observations is the best of our sample ((/.'ST). The effect of this 
excellent PSF is most obvious in the central velocity field of the 
ionised gas. The large-scale stream of gas extending to the centre 
from larger radii is clearly decoupled from a rapidly rotating central 
disk component. This disk of gas is well known from narrow band 
imaging studies, also from HST fe.g. lFord et alll994l) . A closer ex- 
amination of Fig. ll2l shows that the velocity dispersion of the gas 
in this galaxy is also very high in the central parts, suggesting the 
unresolved superposition of two or more velocity components. Van 
der Marel 1 1994) also found a large central velocity dispersion of 
the H7 emission line, although his value of 516 km s~^ is signifi- 
cantly lower than our central value based on [OIII] (> 700 km s^ ^ ). 
This difference can be attributed to the higher spatial resolution of 
our data (0'.'57 FWHM sampled with 0'.'27 lenses, compared with 
0'.'79 FWHM sampled with a 1" long-slit). Fig. [T3lp resents the cen- 
tral spectrum of the OAS I S data, showing that the observed emis- 
sion lines are mainly composed of two strong components, belong- 
ing to the approaching and receding sides of an unresolved disk, 
where the absolute velocity amplitude is over 1000 kms"\ This 
is in agreement w ith spectroscopic measurements with HST (e.g. 
iHarmset all 1994 . and is thought to indicate emission from a gas 
disk orbiting in the Keplerian potential of a 3 x 10^ Mq black hole 
(Macchetto et al. 1997). 

Another interesting case is the central disk in NGC 3379. This 
disk appears to be rather misaligned with the stellar kinematics, 
although is very well aligned with the cent ral dust ring visible i n 
the unsharp-masked HST image (see also iGebhardt et ai]|2OO0t) . 
The nature of this central gas structure is explored in detail by 



IShapiro et all j200d) . combining our OASIS data with ST IS spec- 
troscopy. 

The gas disks in our sample are generally quite faint and 
difficult to detect, as well as being restricted to the very central 
regions. By contrast, several galaxies show much higher emission- 
line flux distributed across the whole OASIS field. The velocity 
field morphology of these objects generally shows a twisted 
structure, where the axis of rotation smoothly changes with radius. 
Galaxies showing this behaviour include: NGC 2768, NGC 2974, 
NGC 3032, NGC 3414, NGC 3489, NGC 4150, NGC 4486 (in the 
outer regions), and NGC 5846. Several of these objects also show 
an elongated structure in the gas, which tends to be almost aligned 
with the rotation of the object, indicating either a close to edge-on 
disk (from the high observed flattening) or perhaps a bar structure. 

Galaxies exhibiting clearly detectable emission, but which 
have gas that is irregularly distributed or does not show a clear 
sense of rotation are NGC 2549, NGC 4262 and NGC 5982. In 
NGC 2549 and NGC 5982, the gas structure indicates the presence 
of streams of material. This explanation is supported by the 
SAURON observations, which show these structures extending to 
larger radii. In the case of NGC 4262, the physical explanation for 
this lopsided structure is less clear. The SAURON data show ionised 
gas across most of the field of view, indicating a significant amount 
of gas, rather than a relatively small filament. The central regions, 
however, show a large asymmetric feature, which connects to the 
same feature in the OAS I S data. The unsharp-masked HST image 
also shows a peculiar ring-like feature, with the signature of faint 
dust lanes across it, co-spatial with the strongest emission features. 
This galaxy exhibits a strong stellar bar, the dynamical influence of 
which may help explain these peculiar features. 

Three objects have detectable emission lines, but have a very 
'patchy' distribution, generally exhibiting emission close to our 
detection limit: NGC 1023, NGC 4459, NGC 4621. These galaxies 
also show evidence for emission in the SAURON data, albeit at lower 
spatial resolution. Due to the patchy distribution in our OASIS 
data, the rotation structure in these objects is difficult to determine, 
although in the case of NGC 4459, rotation follows the stars. 

4.2.3 Ionised Gas Kinematics: Dispersion 

Even with our detailed technique for separating the emission and 
absorption components of our spectra, measuring the velocity dis- 
persion of the gas is much more sensitive to noise and template mis- 
match than the flux or velocity. Thus, in the case of weak emission, 
there is a large scatter in the determined dispersion values, with a 
tendency to measure broad lines (Paper V). A more difficult prob- 
lem occurs with galaxies that have large stellar velocity dispersions, 
in that small discrepancies between our template library and the ac- 
tual galaxy stellar population caused by differences in abundance 
ratios give rise to template mismatch effects, which can be com- 
pensated by including an artificially broad emission component in 
the fit. Where emission is strong enough, however, the width of the 
emission lines can be measured reliably without strong biases. 

In general, the emission-line dispersion maps show values be- 
low that of the stars, but with a strong central peak, which can be 
as large or larger than the central stellar dispersion at the same po- 
sition (e.g. NGC 2974, NGC 3414). For most objects, this central 
peak represents unresolved rotation around the galaxy nucleus. This 
effect is most clear in NGC 4486, where the measured velocity dis- 
persion becomes enormous (> 700 km s~^). As mentioned above, 
this comes from fitting two lines, which are separated by around 
1000 km s^^, with a single Gaussian. This is an extreme case: no 
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other examples of clearly separate lines are found, and we consider 
a single Gaussian to be a good representation of the emission lines 
in our data. 

Where the gas rotation is regular and spatially well-resolved, 
the dispersion is generally low (e.g. NGC2768, NGC3032, 
NGC3414, NGC3489, NGC4459, NGC4526 and NGC5198). In 
this respect, NGC 2974 is peculiar, in that the rotation field appears 
rather regular, but the disper sion is q uite high within most of the 
OAS I S field. At lai-ger radii, iKrainovic et alJ <2005l) were able to 
model the gas in this galaxy with near-circular motion under the 
assumptions of asymmetric drift, treating the gas as collisionless 
'clumps' moving under the influence of gravity. In the central re- 
gions, this treatment was still unable to reproduce the spatially- 
extended high velocity dis persion of the gas there. From TIGER 
integral-field spectroscopv. Emsellem. Goudfrooii & FerruilH2003h 
also found an increase in the central emission line widths, and ac- 
count for this as non-circular motions giving rise to a complex line 
profile composed of multiple components, indicating the presence 
of a bar structure. Fitting this with a single Gaussian gives rise to 
the high broadening we observe. This structure is visible in their Ha 
spectral region emission lines (spectral resolution ~ 3 A), although 
we do not detect statistically significant departures from Gaussian 
profiles in the fainter blue emission lines in our lower-spectral res- 
olution data. 



4.2.4 [0III]/HP Ratios 

The two strongest emission lines in the wavelength region used, 
namely H/3 and the [OUT] doublet, may originate from different 
kinds of gas clouds along our line-of-sight, excited by different ioni- 
sation mechanisms, and emitting different amounts of H/3 and [OIII] 
photons. The ratio of the fluxes of these two lines [OlII]/H/3 holds 
information about the physical state of the gas, and in particular 
on its degree of ionisation. These two lines alone are not sufficient, 
however, to understand the role of the various ionising mechanisms, 
although very low values of [OlII]/H/3 (< 0.33) are generally ob- 
served only in the presence of star formation. Indeed, H/3 emission 
generally dominates over [OIII] in star-forming regions, unless the 
metallicity of the gas is sufficiently low. As in Paper V, we consider 
it unlikely that the metallicity of the gas varies abruptly within a 
given galaxy, although it may vary between objects, should the gas 
have a different origin or enrichment history. 

The following discussion is restricted to the 14 objects exhibit- 
ing clear evidence of both emission lines. Within this subset, there 
are a large variety of global ratio values and local structures. From 
the maps in Fig. 1121 the clearest connection between the [OlII]/H/3 
ratio and other galaxy properties is with the dust. Similarly to 
Paper V, we find that the galaxies with the lowest [OlII]/H/3 ratio 
(NGC 3032, NGC 4459 and NGC 4526, which afl have [Om]/H/3 
significantly less than unity) also show extensive dust in the form of 
a regular and settled disk. The other objects which show clear signs 
of dust extinction in their unsharp-masked HST image (NGC 2768, 
NGC 2974, NGC 3379, NGC 3489, NGC 4150, NGC 4262, 
NGC 4486 and NGC 5813), however, have [OlII]/H/3 values around 
unity, and in some cases, significantly larger In these latter objects, 
the dust is either less significant, or has an irregular (but apparently 
planar) distribution. This suggests that star-formation occurs more 
readily where the dust indicates cold, relaxed disk-like structures. 

Although we do not have sufficient emission-line diagnostics 
to interpret our line ratios directly, in some cases there is additional 
evidence that very low [OlII]/H/3 values indicate star-formation. In 
NGC 4526, the regions of strong dust extinction closely correspond 



to regions of low [OlII]/H/3 values, which in turn correspond to re- 
gions of young stellar populations. The metallicity of the stars in 
these regions also appears to decrease, indicating the gas is proba- 
bly not very metal-rich, thus supporting our interpretation that star- 
formation is the dominant ionisation mechanism in those regions. 
The counter-example to this is NGC 3489, which also shows dust 
and young stellar populations of moderate metallicity. This galaxy, 
however, exhibits the highest [OlII]/H/3 in our sample, illustrating 
the limitations of interpreting [OlII]/H/3 on its own. Clearly addi- 
tional diagnostics, such as Ha and [Nil], would be valuable in dis- 
criminating between the true ionisation mechanisms in our galax- 
ies, and establishing the role of dust and gas dynamics in the star- 
formation process. 

4.3 Line Strengths 

Figure fT2l includes maps of the six Lick indices within our OAS I S 
spectral range. We measure four separate iron indices, of which 
Fe5015 is affected by [OIII] emission. It is reassuring that, even 
in galaxies with significant emission, the morphology of the differ- 
ent iron index maps is generally consistent (e.g. the central peak in 
NGC 3032, or the dusty region in NGC 4526). The two other mea- 
sured indices, namely H/3 and Mg b, are often used as diagnostics 
of stellar population age and metallicity (respectively), usually in 
combination with additional iron features. 

Applying the method described in Section lSiSl we also present 
maps of luminosity-weighted stella r age, metallicity and abun- 
dance ratios using the SSP models of|Th omas. Maraston & Bended 
i2003h . Fig. I14| presents the age, metallicity [Z/H] and abundance 
ratio [a/Fe] derived from the 4" aperture spectrum described in pre- 
vious sections, each plotted as a function of the stellar velocity dis- 
persion measured inside the same aperture. The values and errors 
(marginalised over the other two free parameters) of these parame- 
ters are given in TableQ as is the fraction of the galaxies' effective 
radii covered by the 4" aperture. Given our limited field of view 
and the large range in effective radius of our targets, no attempt was 
made to account for the different intrinsic aperture sizes. We per- 
form a robust straight-line fit to the data, the parameters of which 



log* = -3.95 + 2.04 log cr, 
[Z/H] =0.02 + 0.14 log a, 
[a/Fe] = -0.45 + 0.26 log a. 



(3) 
(4) 
(5) 



As found by many authors (e.g. |Ca rollo^Danzi£er. & Buso j 
1993 jTrager et a lj2OO0l:rKuntschner et alJ200lllProctor & Sansoij 
2002l: Bernardi et alj |2005t iThomas et alJl2005l> . we find that the 



stellar population parameters correlate with velocity dispersion. 
This is in addition to the empirical trends of individual indices 
and combinations of indices with velocity dispersion, such as 
the well-known tight correlation with Mg absorption strength 
(e.g. Bernardi et al. 1998; Colless et al. 1999; Worthev & Collober^ 
I2OO3L and others). We find rather a strong trend with age, such that 
lower-dispersion galaxies show younger luminosity-weighted ages. 
This is consistent with th e 'downsizing' scenario of galaxy evo- 
lution iCowie et alJll99^ . where the characteristic mass of star- 
forming galaxies decreases with time. Our small magnitude-limited 
sample is biased by some degree at lower velocity dispersions 
(masses) towards young objects, resulting in a steeper trend. Given 
the small obser ved range in mass-to- light ratio (M/L) with velocity 
dispersion (e.g. lBemardi et a l.''2003''). however, this effect is small, 
and will not change the trend of decreasing age with dispersion. We 
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do note, however, that w e derive a significantly stronger gradient 
than lThomas et alj <2005h . 

A correlation also exists with [a/Fe], such that higher- 
dispersion galaxies have larger [a/Fe] values. This trend has been 
inferred empirically from index ratios (e.g. ?), and also measured 
directly (as here) using population model s with variable abun - 
dance ratios (e.g. *Proctor & Sanso 1^ 12001 Irto ^as et alJl200.it) . 
This trend implies a variation in star- formation timescales, such that 
higher-dispersion (higher-mass) systems have experienced more 
rapid star formation than lower-dispersion galaxies. Linking to the 
trend with age, this rapid star-formation phase happened at early 
epochs, leading to the existence of massive star-forming systems 
at high redshift, perhaps explai ning th e existence of a population 
of bright sub-mm sources (e.g. fSmail e t al. 2002). We find only a 
weak relation between decreasing metallicity and velocity disper- 
sion, and there is a broad range of metalli city values for any given 
velocity dispersion. Other authors (e.g. lThomas et al]|2005l) find a 
clearer relation, albeit in the same sense as ours. 

The inhomogeneous aperture and small sample size strongly 
limit the generality of trends we find from these aperture measure- 
ments, although they are in broad agreement with other authors. In 
the following sections, we focus on the structure found in the two- 
dimensional maps of population parameters. 



4.3.1 Luminosity-weighted age 

Most objects in our sample show almost no significant variations 
in age within the OASIS field of view. There are a number of 
notable exceptions to this, however, which show sharp changes in 
the spatial distribution of stellar ages. Three galaxies in our sample 
show a clear decrease in age towards the galaxy centre (NGC 3032, 
NGC3489 and NGC 4150). The most extreme case is NGC 3032, 
which reaches a minimum age of 0.5 Gyr in the central OAS I S bin. 
These sharp drops in age derived from the models are driven by 
sharp increases in the H/3 line-strength clearly visible in the maps. 
Also, in both NGC 3032 and NGC 3489, the peak in H/3 absorption 
is coincident with a decrease in the H/3 emission equivalent width, 
indicating that the observed strong H/3 absorption is not an artifact 
of overestimated emission subtraction. 

Another interesting case is NGC 4526. The age map clearly 
shows a young population associated with the dusty disk embedded 
in the central 10" of this galaxy. The foreground (northern) side of 
the disk shows the age decrease most clearly, although the struc- 
ture can be traced also on the other side of the bulge. There is a 
sharp rise in stellar age at the northern edge of the field, which can 
be associated with the corresponding edge of the dust-lane in the 
unsharp-masked HST image. 



4.3.2 Metallicity 

All galaxies in our sample show constant or increasing metallic- 
ity towards the galaxy centre, although the variations are generally 
small within the OASIS field. In Paper VI, we found that flattened, 
fast-rotating galaxies show an Mg b distribution which is flatter than 
the apparent surface brightness in the SAURON data, suggesting the 
presence of a metallicity- or abundance ratio-enhanced disk em- 
bedded in these objects. We do not expect to find strong evidence 
of this effect in our OASIS maps, since the disk component is not 
expected to dominate over the spheroidal component in the cen- 
tral regions. For some of the fast-rotating, flattened objects, how- 
ever, we do find weak evidence of a metallicity distribution that is 
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Figure 14. Relations between stellai' velocity dispersion a and age, 
metallicity and abundance of our galaxy sample, measured using the 
Thomas, Maraston & Bend^ stellar population models. Values and 

errors are given in Table Q The models are fitted to the six Lick indices 
in our spectral range, measured from the 4" radius circulai' aperture inte- 
grated spectrum described in previous sections. Symbols indicate the RC3 
classification: circles = ellipticals, triangles = lenticulars. Filled symbols in- 
dicate the cluster environment sample, open symbols indicate field galaxies. 
Oveiplotted are robust straight line fits to our data, with parameters given in 
eQS. I3l5l . Dotted horizontal lines indicate the limits of the SSP model range. 



flatter than the surface brightness distribution (namely, NGC 3384, 
NGC 4473 and NGC 5845). Of these, NGC 5845 is pai'ticularly in- 
teresting. The stellar kinematics indicate the presence of an embed- 
ded disk (Section l4.1.2t , which is clearly visible in the unsharp- 
masked HST image. In the map of metallicity, there is evidence for 
a weak enhancement along the major axis within the central 2", 
which may be associated with this thin disk component. 

In the three galaxies with young nuclei mentioned above, the 
young population is associated with an increase in the apparent 
metallicity, implying that the young stars (which are dominating 
the light in these regions) are more metal rich than those at larger 
galactic radii. By contrast, within the young disk of stars associ- 
ated with the dust structure in NGC 4526, the metallicity shows a 
distinct decrease compared to the metal-rich bulge region. The ab- 
solute metallicity of the newly formed stars in all cases is, however, 
similar, and generally super-solar. 
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Table 7. Central stellai' population parameters, determined from six Lick in- 
dices measured on a single OASIS spectrum integrated within a 4" circular 
aperture. 



NGC 
(1) 


Re 

(2) 


a 
(3) 


Age 
(4) 


[Z/Fe] 
(5) 


[o/Fe] 
(6) 


1023 


0.08 


210 ±6 


4 7+3.6 
^■'-2.0 


0.52 ±0.12 


0.14 ±0.05 


2549 


0.20 


142 ±5 




0.37 ±0.11 


0.09 ± 0.04 


2695 


0.19 


225 ±4 




0.22 ± 0.08 


0.09 ± 0.05 


2699 


0.29 


132 ±7 


^■•'-0.2 


0.67 ± 0.04 


0.21 ± 0.05 


2768 


0.06 


191 ± 3 




0.60 ±0.11 


0.16 ±0.04 


2974 


0.17 


236 ±5 


7 9+3.7 
'■^-2.5 


0.15 ± 0.04 


0.07 ± 0.05 


3032 


0.24 


107 ± 8 


9+0-5 


-0.08 ± 0.38 


-0.13 ±0.11 


3379 


0.10 


224 ±6 


11 5+^-^ 


0.22 ±0.12 


0.18 ± 0.08 


3384 


0.15 


162 ±5 


3 2+1-5 


0.44 ±0.11 


0.16 ±0.04 


3414 


0.12 


237 ± 5 


3g+3.7 


0.44 ± 0.22 


0.09 ± 0.05 


3489 


0.21 


98 ±4 


1 7+0.2 
^■'-0.2 


0.22 ± 0.08 


0.03 ± 0.04 


3608 


0.10 


190 ±5 




0.22 ±0.12 


0.14 ±0.04 


4150 


0.27 


91 ±6 




0.07 ± 0.08 


0.07 ± 0.06 


4262 


0.40 


174 ±5 


10.21^:^ 


0.22 ± 0.08 


0.12 ±0.05 


4382 


0.06 


179 ±6 


-1 7+0.3 
^■'-0.2 


0.52 ±0.15 


0.12 ±0.06 


4459 


0.11 


169 ±5 


1 9+0.1 


0.67 ± 0.04 


0.14 ±0.04 


4473 


0.15 


192 ±4 


13 2+*-° 


0.07 ± 0.12 


0.16 ±0.05 


4486 


0.04 


361 ± 9 


15 1+^-1 


0.22 ± 0.04 


0.19 ±0.03 


4526 


0.10 


247 ±6 


7 0+0.2 


0.67 ± 0.04 


0.25 ± 0.05 


4552 


0.13 


281 ± 5 




0.44 ±0.11 


0.23 ± 0.05 


4564 


0.19 


156 ±4 


3 f;+3.3 
^■"-1.7 


0.52 ±0.15 


0.12 ±0.05 


4621 


0.09 


233 ± 6 




0.52 ± 0.08 


0.16 ±0.05 


5198 


0.16 


201 ± 8 




-0.08 ± 0.04 


0.30 ± 0.04 


5813 


0.08 


240 ±4 


15 1+^1 


0.15 ± 0.04 


0.19 ±0.03 


5831 


0.11 


158 ±4 


T, 7+3.4 
^■^-1.7 


0.44 ± 0.22 


0.14 ±0.06 


5845 


0.87 


237 ± 5 


gQ+2.6 


0.30 ±0.11 


0.09 ± 0.07 


5846 


0.05 


247 ± 7 


7 9+5-2 
'■^-3.2 


0.44 ±0.11 


0.16 ±0.05 


5982 


0.15 


257 ±6 


1 7+2.8 
i-'-^_2.3 


0.00 ± 0.08 


0.12 ±0.04 



Notes: (1) NGC number (2) Fraction of effective radius covered by 4" 
aperture. (3) Velocity dispersion measured within 4" aperture, in km s~^. 
(4) Luminosity-weighted age measured within 4" aperture, in Gyr. (5) [Z/H] 
measured within 4" aperture. (6) [cj/Fe] measured within 4" aperture. Er- 
rors are given as (one-dimensional) Icr, marginalised over the other two 
parameters. 

4.3.3 Abundance ratio 

With the exception of two galaxies, the maps of abundance ra- 
tio either show a weak increase towards the centre, or are rather 
flat, showing little significant structure within our field of view. 
NGC 4150 shows a weak decrease in [a/Fe] towards the centre. 
In this case, the decrease is also coincident with younger stars and 
a decoupled kinematic component. This is also true for the disk of 
young stars in NGC4526, where the abundance decreases within the 
dusty region. With the limitations of interpreting integrated popula- 
tions with SSP models in mind (see l3.8> . given the similar spatial 
structure in the maps of age and abundance for these two galax- 
ies, it seems that the decrease in apparent abundance is due to the 
presence of the young stars. 



NGC 3032 shows a particularly low [a/Fe] in the outer parts of 
our field, giving the only sub-solar abundance ratio in our sample. 
In the central regions, however, the abundance becomes near-solar 
coincident with the strong decrease in age. This relative enhance- 
ment of ci-elements with respect to iron suggests that the young 
stellar population is strongly enriched with material expelled from 
short-lived massive stars via type II supemovae. This implies that 
either the star-formation timescale within the central component is 
markedly shorter than in the main body of the galaxy, or that the ini- 
tial mass function (IMF) is skewed towards higher-mass star in this 
regi on. A combinatio n of these effects is also a possible explanation 
(e.g. lMatteuccJl994 . 

Several authors have reported that KDCs show enhanced 
metallicity and/or abundance, inferred from changes in line- 



strength gradients coincident with the kinematic structure (e.g. 


ISender & Surma 1992: ICaroUo & Danzieer 


1994b; 


Morelli et alj 



|2004). We note that two of the objects mentioned above, NGC 4150 
and NGC 3032, contain stellar KDCs, but show contrasting changes 
in abundance ratio within the KDC. These KDCs also show dis- 
tinctly young ages (see following Section), unlike the objects stud- 
ied previously, and so may not be comparable. Full consideration of 
this phenomenon for our sample requires also the use of SAURON, 
since most of the KDCs are larger than the OAS I S field. We there- 
fore leave further discussion of this issue to future papers of this 
series. 



5 DISCUSSION 

5.1 Stellar kinematics of young galaxy centers 

From Paper VI, SAURON observations show that galaxies exhibit- 
ing the strongest global H/3 absorption (measured within one effec- 
tive radius) also tend to show an increased H/3 line strength towards 
the central regions. Three of the four galaxies with global H/3 ab- 
sorption around 3 A or stronger are in our OAS I S sample, namely 
NGC 3032, NGC 4150 and NGC 3489. As discussed above, these 
objects show a sharp decrease in their central luminosity-weighted 
age in our OAS I S field, indicating very young stellar components 
that are in some cases not spatially well-resolved, even at the spatial 
resolution of our data. Although these galaxies also show young or 
intermediate aged stars across their main body, their central regions 
(< few 100 pc) are significantly younger. These young centers are 
responsible for driving the spread in age covered by the SAURON 
sample galaxies to include relatively low ages (< 2 Gyr). 

From HST imaging studies, a number of early-type galax- 
ies are found to contain point-like nuclei, unresolved even at 
//irresolutions (^5 pc; e.g. Carollo et al. 1997a; Rest et al. 2001; 
iLaine et al."2003V Lauer et al . 2005; Cote et al. 2006). These nuclei 
are generally blue in broad-band colours, and have been attributed 
to both non-thermal sources and nuclear star-clusters. We differ- 
entiate these nuclei from the young centers found in our data, as 
we do not have sufficient spatial resolution to directly connect the 
H/3-enhanced structures to such small components. Also, the nuclei 
found with HST exist in early-type galaxies of a range of masses, 
whereas the young centers found in our data are confined to the 
low-mass objects. A second important difference is that the central 
structures found with our data are found to be young spectroscopi- 
cally, thus largely avoiding the problems of the age-metallicity de- 
generacy and dust absorption that often hamper the interpretation 
of broadband colour gradients. We do not exclude, however, that 
the young central structures we observe may correspond to nuclei 
at 5 pc scales. 
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As well as the unresolved blue nuclei found with HST, more 
extensive blue components are often seen in so-called 'E+A' galax- 
ies: objects which are morphologically similar to early-type galax- 
ies, but exhibit s trong Bal mer absorp tion typical of A-stars (e.g. 
iDressler & Gunnllll983: Zabludoff et al.^1996) . Moreover, compact 
blue comp onents are found i n these systems, unresolved at scales of 
~ 100 pc I'Yang et alJl2004) . These are thought to be star clusters 
formed during a recent star-burst. Their locations within the main 
galaxy can vary, however, and they are not always confined only to 
the central parts. E-l-A galaxies are a somewhat mixed class of ob- 
jects, whose members often show disturbed morphologies at large 
scales, suggesting a recent major merger The morphologies of our 
young-center galaxies do not show obvious signs of a recent major 
merging event (e.g. tails, shells). However, the compact nature of 
our young components in our galaxies suggests that they may have 
experienced a similar mode of localised star-formation. 



Clusters of recently formed stars are often found 
centers of later type than E/SO (e.g. spirals and late-t^ 
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11998.') . The best candidate in our sample for such a cluster is 
the youngest object, NGC3032, an SABO galaxy that shows a 
bright central peak and possible circum-nuclear ring structure in the 
unsharp-masked HST image (see Fig. I12t . which could be form- 
ing stars. For the next-youngest objects, NGC 3489 (SABO) and 
NGC4150 (SO), this explanation is less compelling for their young 
centers, given the lack of bright central sources or star-forming 
rings, but this may be due to projection and dust. In any case, both 
objects exhibit recent circum-n uclear star-formation, and are clas- 
sified as 'transition' objects bv lHo. Filippenko. & Sar gent ( 1997), 
having emission-line properties betw een that those of HII nuclei 
and LINERS (Ho, FilioDenko, & Sarg enll993l) . 

It has long been suggested that some early-type galaxies show 
young apparent ages du e to a small fraction (by ma ss) of young 
stars within the galaxy <Worthev et ai]ll994 iTraser e t al. 2000). 
This 'frosting' of young stars has a low mass-to-light ratio; thus 
while the young stars do not contribute much in terms of mass, 
their contribution to the integrated light is significant. However, it 
was not clear where these stars were located in the galaxy, or what 
triggered their recent formation. In Paper V we have shown that 
many of our sample galaxies contain large quantities of (ionised) 
gas that may exist in kiloparsec scale structures, suggesting that 
there are extended reservoirs of star-forming material available in 
these galaxies. It would seem, however, that only in the central re- 
gions, deep in the galaxy potential, are the conditions suitable to 
trigger efficient conversion of this material into stars. 

Investigating the stellar kinematics of the three youngest 
galaxies in our sample, the two with the strongest H/3 absorption 
(NGC 3032, NGC 4150) harbour small KDCs, which are essentially 
aligned but counter-rotating with respect to the main body of the 
galaxy (APAkiii~ 180° and ^ 190° respectively). The other galaxy 
(NGC 3489) does not show evidence for counter-rotation, but does 
exhibit a mild twist in its velocity field (APAkin~ 10°) within 
a 2" radius, as discussed in Section l4m Furthermore, the gas in 
this galaxy shows an elongated structure with enhancements at ei- 
ther end, and open spiral features in the dust extinction (Fig. I12> . 
suggesting the possible presence of a central bar (cf. NGC 2974, 
lEmsellem. Goudfrooii & Ferruitl2o63l) . 

In all three galaxies, the region showing structure in the veloc- 
ity field is coincident with the young stellar populations, implying 



that the two phenomena are directly linked. Fig. llSl illustrates this, 
showing the radial profiles of both the kinemetric analysis (see Sec- 
tion and the stellar population parameters. The clearest asso- 
ciation between the stellar kinematic structure and the population 
parameters is for NGC 3032. The 180° jump of the kinematic PA at 
a radius of ~ 2" and the corresponding minimum in rotation be- 
tween the outer body and the central KDC coincide closely with a 
downturn in the age, and upturns in the metallicity and abundance 
ratio. NGC 4150 shows less dramatic variations in the kinematics, 
but is strongly affected by poor seeing. However, PAkin changes 
significantly within the central few arcseconds, and the mean rota- 
tion is almost zero. This region corresponds to the KDC, and also 
shows a decrease in age compared to the outer regions. There is 
also weak evidence of an increase in metallicity and a decrease in 
[a/Fe]. NGC 3489 does not exhibit a counter- rotating component. 
Fig. ^1 shows, however, that PAkin changes by ~ 10° within the 
region where the age decreases from 2 to 1 Gyr The metallicity 
also increases more steeply within this region. In all three cases, a 
change in the stellar populations, and most notably a decrease in the 
luminosity-weighted age, is accompanied by a change in the stellar 
kinematics, suggesting a link between the two. 

5.2 Stellar ages of KDCs 

In the above discussion, we highlighted that two of the three 
youngest galaxies in our sample have a KDC associated with 
the young stars. Here we extend this discussion to examine the 
stellar populations of all KDCs that have been found so far in the 
SAURON survey, considering the combined information from both 
the SAURON and OASIS data wherever possible. 

With integral-field spectroscopy, we can obtain a complete 
two-dimensional map of the KDC, allowing us to measure accu- 
rately the orientation and spatial extent of the structure, as well 
as determining the stellar populations. The SAURON representative 
sample contains numerous examples of KDCs (Paper III), many of 
which are also in our present subsample. To these, we add the three 
KDCs which are visible in our OASIS data, though not resolved 
with SAURON: NGC 3032, NGC 4382, and NGC 4621. 

Table m lists the galaxies from the full SAURON sample that 
we consider to harbour a KDC. These objects are selected as they 
exhibit subcomponents that either rotate around a different axis to 
the main galaxy, or that rotate around the same axis, but with an op- 
posite sense of rotation. With this definition, we do not consider 
dynamically decoupled structures which co-rotate with the main 
galaxy, such as nuclear disks, even though there are several such 
structures in our sample. Misaligned and counter-rotating compo- 
nents offer the advantage, however, of being clearly distinguished 
from the rest of the galaxy, and may indeed have a distinct forma- 
tion process from the co-rotating case. 

Considering both the OASIS and SAURON data, there is a 
large range of apparent sizes for the KDCs in the SAURON sam- 
ple. We estimate the size of the KDC by using the isolating region 
in the velocity map which surrounds the KDC, where the superpo- 
sition of the two components gives a local minimum in the mean 
velocity. Although approximate, this gives a straightforward esti- 
mate of the projected extent of the central structure. We also include 
three galaxies which show evidence of a KDC in the SAURON ve- 
locity fields of Paper III, but for which no OAS I S observations were 
obtained (NGC 4458, NGC 7332 and NGC 7457). In addition, we 
include the well-known KDC in NGC 4365, which was observed 



with SAURON, though not as part of the main survey I Davie s et alj 
l200ll) . We re-measure the KDC diameter for this galaxy using our 
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Figure 15. Radial profiles of relative kinematic position angle (APA]ji,j); maximum rotation velocity, as given by the fci term of the kinemetric analysis; and 
stellar population parameters: age, metallicity (Z) and abundance ratio ([o/Fe]). Small dots indicate the individual stellar population parameters measured from 
our index maps. Red diamonds indicate the biweight mean value of these individual measurements within a circular annulus, plotted against the geometric 
radius of this annulus. The galaxy name is indicated in the title of each panel. 



metho d, and take the age and H/3 hne-strength from lPavies et alJ 
jZOOll) . For each KDC, the data with the smallest PSF and with 
a large enough field of view to encompass the isolating ZVC was 
used to estimate the diameter. The apparent diameters are converted 
to intrinsic sizes using the surface brightn ess fluctuation (SBF) dis- 
tance measurements of lTonrv et aljJlOOll) . Distances, and apparent 
and intrinsic KDC diameters are listed in Table|8| 

Table|8|also lists the central value of the H/3 absorption-line in- 
dex measured from a spectrum integrated within a 1" circular aper- 
ture centred on the galaxy, using OAS I S data where possible. This 
index was combined with Mg b and available iron indices within the 
same aperture to estimate the luminosity-weighted age, as in Sec- 
tion]^ For SAURON age estimates, only the Fe5015 and Fe5270 
iron indices are available. 

The values given in Table |8| cover a very broad range of pa- 
rameter space, spanning more than an order of magnitude in both 
size and luminosity-weighted age. Fie. ll6l illustrates this, showing 
that the large 1 kpc) KDCs have predominantly old 10 Gyr) 
populations. The group of small (< 1 kpc) KDCs, by contrast, ex- 
hibit populations spanning the full range of ages, from < 1 Gyr to 
15 Gyr, with five of the six small KDCs younger than 5 Gyr. The 
colours used in this figure indicate the class of rotation (also given 
in Table|8}, being either a fast-rotating object, mainly supported by 
rotation (blue symbols); or a slow-rotating object, mainly supported 
by dynamical pressure (red symbols). There are clearly no examples 
of slow-rotating objects with young, compact KDCs. 

In addition, the small, young KDCs observed with OASIS 
are all very close (i.e. within 10-20°) to being counter-rotating 
(i.e. APAkiii~ 180°), as far as our data can reveal. In the cases 
of NGC4150 and NGC4621, our measurements of PAkin for 
the KDC are strongly affected by seeing effects. For NGC3032 
and NGC4382, however, the measurements are robust. The other 



two young KDCs (NGC 7332 and NGC 7457) are only observed 
with SAURON, and limited spatial information prevents a detailed 
measurement of the KDC PAkin- By subtracting a model of the 
large-scale rotation, Falcon-Barroso et al. (2004) estimate the 
PAkin of the KDC in NGC 7332 as 155° (25° from the photometric 
major-axis), consistent with our general finding of near-aligned 
counter-rotation. The small misalignment of these KDCs raises 
the question of the existence of co-rotating counterparts. Such 
components have a more subtle effect on the velocity field, having 
similar rotation amplitude to the background large-scale rotation. 
For this reason, we limit the current analysis to the clearly visible 
KDCs, and note that the majority of young nuclei in our sample are 
in fact KDCs of this kind. 

An indicator of how significant the KDC is in terms of the 
global galaxy properties is the fraction of the galaxy's mass which 
is contained in the KDC. Accurately estimating this quantity is dif- 
ficult. Perhaps the most accurate but complex way is via dynamical 
modelling, where the KDC can be separated in phase-space directly 
(e.g. Cap pellar i et al. 2002). Such work is beyond the scope of this 
paper, but will be addressed in a future paper of this series. A very 
simple way to quantify the significance of the KDC is to consider 
its size with respect to the effective radius Re. The young KDCs 
generally extend to only < 0.1 i?e, whereas the older KDCs are 
generally 0.2 - 0.3 i?e in size. Moreover, the young KDCs have a 
low mass-to-light ratio, thus contributing significantly to the light, 
but relatively little in mass. 



5.3 Age distribution of compact KDCs 

Fig. 1161 shows that the compact (i.e. smaller than 1 kpc diameter) 
counter-rotating KDCs occur only in the fast-rotating galaxies of 
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Table 8. List of galaxies with KDCs from SAURON and OAS I S stellar kine- 
matics. 



NAME 


Dist. 


Diameter 


Diameter 


Age 


Rotation 


(NGC) 


(Mpc) 


(") 


(pc) 


(Gyr) 




3032 


22.0 (3.0) 


3.5 


372 (73) 


0.50 (0.1) 


F 


3414 


25.2(4.1) 


23.0 


2813 (465) 


18.00(2.3) 


S 


3608 


22.9(1.5) 


21.0 


2332 (164) 


17.00(1.0) 


S 


4150 


13.7(1.6) 


4.0 


266 (45) 


1.00 (0.1) 


F 


4382 


18.5 (1.2) 


3.5 


313 (49) 


3.00 (0.3) 


F 


4458 


17.2(1.0) 


10.0 


834 (63) 


17.00(1.5) 


S 


4621 


18.3 (1.8) 


2.0 


17/ (47) 


1 O f\l\ /"I 1 \ 

18.00 (2.1) 


F 


5198 


36.3 (3.5) 


6.0 


1056(134) 


15.00 (2.5) 


S 


5813 


32.2(2.8) 


20.0 


3123 (280) 


17.00 (2.5) 


S 


5831 


27.2 (2.2) 


18.0 


2370 (203) 


9.00(1.2) 


S 


5982 


41.9(4.0) 


10.0 


2030 (220) 


8.00(1.5) 


S 


7332 


23.0 (2.2) 


4.0 


446 (70) 


3.50 (0.2) 


F 


7457 


13.2(1.3) 


4.0 


256 (41) 


2.00 (0.2) 


F 


4365 


20.4(1.7) 


15.0 


1484 (130) 


14.00 (2.5) 


S 



Notes: Each column gives (from left t o right): galaxy name; distance from 
surface brightness fluctuations of Tonry et al. (2001); estimated apparent 
diameter from SAURON or OASIS velocity field; intrinsic diameter in pc; 
luminosity-weighted age estimate from hne-strengths measured withi n 
1" using the population models of iThomas. Maraston & BendeJ |2223)i 
galaxy classification according to the measured value of specific angular 
momentum within one effective radius (Paper IX, in preparation). Errors 
for each parameter are given in parentheses. Errors in distance come from 
and the eiTor in diameter combines the distance en'or 
with an estimated accuracy of 0^'5 for measuring the KDC diameter from 
the velocity map. 



our sample. Figure fTTl shows that the distribution of fast-rotating 
galaxies with counter-rotation in our sample is skewed towards ob- 
jects with younger global luminosity-weighted mean ages. Five of 
the seven objects younger than < 5 Gyr show a counter-rotating 
core (and one of these has still to be observed with high spatial res- 
olution, and so may yet reveal a KDC). Only one of the remaining 
thirteen objects older than 5 Gyr shows a detectable KDC. 

The higher occurrence of counter-rotation in young galaxies 
may simply be the effect of luminosity weighting. These small com- 
ponents contribute relatively little in mass to the galaxy, but are still 
visible above the background of the main body because their young 
age (and therefore low M/L) makes them very bright. We explore 
this with the following 'toy' model. Using the measured central 
line-strengths of NGC 4150, we constrain the amount of mass of 
young stars within the central aperture that can be added to a back- 
ground 'base' population (assumed to be that of the outer parts of 
the galaxy). Our SAURON data indicate an age of 5 Gyr for the 
outer body of this galaxy. The central H/3 value is ~ 3.6 A, allowing 
a mass fraction of around 6% of a 0.5 Gyr population to be added. 
Note that this is only the contribution within the 1" line-of-sight, 
not to the entire system. 

We simulate a two-component velocity field using a kineme- 
try Fourier expansion with ad hoc coefficient profiles, consisting of 
a large-scale rotating component, with a compact counter-rotating 
component. Both components have the same velocity dispersion. 
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Figure 16. Intrinsic size of KDCs versus their centi'al luminosity-weighted 
age. Closed symbols show objects observed by both SAURON and OASIS. 
Open symbols show objects only observed with SAURON. Red points indi- 
cate slow-rotators and blue points indicate fast rotators, as defined in Paper 
IX (in preparation). Arrows indicate where the measured diameter of the 
KDC is less than four times the seeing FWHM, implying that the KDC is 
not well-resolved and may be smaller than measured here. 
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Figure 17. Histogram of ages for the fast-rotating galaxies from the 
SAURON survey. From the complete SAURON sample (black line), we in- 
dicate those observed so far with OASIS (blue), and those with a counter- 
rotating KDC (red). Ages are derived from the SAURON lines-strengths of 
Paper VI. 



We a ssign th e two velocity distribution fields to SSP model spec- 
tra of lBruzua l & Chariot ( 2003), assigning the young population to 
the compact counter-rotating component, and combine the spectra 
using the appropriate mass-weighting. Noise is added to the simu- 
lated data cube, which is then spatially binned, and the kinematics 
are extracted using pPXF, consistent with how we treat our obser- 
vational data. The result is a realistic-looking KDC velocity field 
(Fig. llSt . This test demonstrates that our analysis technique (which 
uses a different library of SSP models as kinematic templates) can 
reliably recover the velocity distribution, showing that the observed 
KDCs are not spurious artifacts of our kinematics analysis due to 
the presence of strong population gradients. 
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Figure 18. Model velocity fields for an evolving counter-rotating (top row) and co-rotating (middle row) young subcomponent. The spectrum of the young 
decoupled component (blue) and older main component (red) are shown on the bottom row, indicating there luminosity-weighted contribution to the observed 
spectrum (black). 



We then hold the mass fraction fixed, and 'evolve' the pop- 
ulations in step by combining older SSP models with the same 
age difference. As the KDC population ages, its mass-to-light ra- 
tio increases, resulting in a dimming of the KDC stars. The ef- 
fect is to 'fade' the KDC into the background rotation field, as 
the luminosity-weighted contribution becomes less significant. Af- 
ter 5 Gyr, the KDC is barely visible. This helps explain the apparent 
lack of intermediate and old aged small KDCs. 

Assuming that young co- and counter-rotating components 
trace the angular momentum vector of the 'seed' material which 
has recently entered the system from outside, both types of compo- 
nent are equally likely to form. It is therefore somewhat surprising 
that five of the seven youngest galaxies in our sample show evi- 
dence of counter-rotating components. These are, however, small- 
number statistics, and considering errors in the age determination 
and what maximum age limit one considers, the significance that 
counter-rotating components tend to be young can be diminished. 
It is tempting, however, to conclude that co unter-rotat ion is some- 
how connected to enhanced star-formation. iKapferer e t al. ( 20oJ) 
explore this through combined N-body/hydrodynamic simulations, 
and find that counter-rotating encounters show a slight tendency to 
produce higher star-formation rates than their co-rotating counter- 
parts, although there are a large number of exceptions to this. 

We have focused our discussion on the counter-rotating case, 
mainly because these components leave an obvious signature in the 
observed velocity field. The middle row of Figure fTsl uses the same 
kinematic components and populations as the top row, but in this 
case the subcomponent is co-rotating. The impact of the co-rotating 
component on the total observed velocity field is clearly more sub- 
tle, visible only as an additional 'pinched' central isovelocity con- 
tour. The co-rotating component becomes difficult to separate from 



the background field after only ~ 2 Gyr. For this reason, the in- 
trinsic distribution of co- and counter-rotating central young com- 
ponents could actually be rather similar, but since the co-rotating 
cases can only be clearly identified at very young ages, our sample 
of young galaxies is currently too small to determine this reliably. 



5.4 Interpretation 

The existence of KD Cs has been kn own since the mid 1980's from 
long-slit studies (e.g. lPranx & Illing worth 1988; Bender 1988), and 
integral-field studies have reinforced this (e.g. iDavies et ali li200I. 
Paper III), allowing their full kinematic structure to be mapped. 
Since the discovery of these structures, they have often been used 
as evidence that coUisionless galaxy mergers have played a sig- 
nificant role in forming early-type galaxies, where the tightly- 
bound core of one galaxy becomes embedded in the centre of the 
remn ant (e.g. .Ko rmendv 1984; Hollev-Bockelmann & Richston^ 
l200(i) . This scenario can naturally account for the often misaligned 
orientation of the KDC, as well as the dynamically hot nature of the 
remnant galaxy. 

An alternative explanation for the formation of KDCs is that 
they 'grow' via star-formation inside the host galaxy. To account for 
misalignment of the KDC, the star-forming m aterial is assumed to 
be accreted from outside the host galaxy (e.g. lHemguist & Barnes! 
1991; Weil & Hernauist 1993). 

In both of these scenarios, the KDC should show distinct stel- 
lar populations from the host galaxy. A KDC 'grown' inside the 
main galaxy would be composed of younger stars than its host. 
A KDC formed as the denser, more tightly-bound core of one of 
the progenitor galaxies of a major merger should also exhibit dis- 
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tinct populations due to the different star-formation histories of 
the two galaxies (which gives rise to the necessary density differ- 
ence). In general, however, little evidence has been found to suggest 
that KDCs have partic ularly distinct star-formation hi stories from 
their host g alaxy (e g. iForbes et alJll996l : ICarollo et aL .1997a.bi: 
lOavies et al1l200ll: lBrodie et alj|2005l but see iBender & Surmal 
119921: iBrown et alJl200 0')." This apparent lack in contrast between 
the stellar populations of the KDC and its host galaxy suggest that 
either the assembly of these objects is rather carefully arranged, 
or that they formed at early epochs, and the observable differences 
between the populations have been diluted by time, masking their 
distinct origins. 

From our high spatial resolution observations, we find a num- 
ber of compact KDCs which do show clear evidence of distinct, 
young stellar populations. Extending this connection between the 
stellar kinematics and populations for all KDCs in the SAURON 
survey using both OASIS and SAURON data, we find two types 
of KDC: large, kiloparsec-scale KDCs which have old stellar popu- 
lations, and are found in slowly-rotating systems; and very compact 
KDCs on 100 pc scales, which are generally very young, are found 
in aligned, fast-rotating objects, and are close to counter-rotating 
around the same axis as the outer parts.. 

This intuitively fits with the picture that fast-rotating galax- 
ies (including possible sub-components) tend to form through dis- 
sipative processes, involving a significant amount of gas. There is 
a sizeable supply of star-forming material for these objects, pro- 
vided it can reach a high enough density to induce star-formation. 
Such material will tend to fall to the centre of these objects, los- 
ing angular momentum via hydrodynamical processes. Reaching 
the critical density for star-formation is most likely at the centre of 
the potential, where the material is collected and compressed. This 
naturally results in compact central star formation. How efficiently 
gas can reach this region may be key to determining the degree of 
star-formation that can occur. In this respect, dynamical perturba- 
tions allow more efficient transfer of angular momentum, resulting 
in a more rapid build-up of material at the galaxy centre and con- 
sequent star-formation. In the presence of such a perturbation (e.g. 
a bar), the orientation of the angular momentum of the gas may be 
significantly altered, which may explain the presence of the stellar 
KDCs having formed from this material. The fact that the KDCs are 
rather well-aligned with the host galaxy supports the interpretation 
that dissipational processes play an important role in the formation 
of these components. 

Slowly-rotating galaxies, on the other hand, are thought to 
form through predominantly dissipationless processes cBinnev 
Il978l: iDavies et alJll983t iBender. Burstein & FabeJll992h . result- 
ing in near-spheroidal systems, often with some degree of triaxi- 
ality. Since these objects have little star-forming material, they are 
less likely to contain young stellar populations. In addition, in a dy- 
namically hot coUisionless system, kinematic subcomponents (ei- 
ther acquired via interactions, or 'grown' from accreted material at 
early epochs) can be dynamically stable over the lifetime of the ob- 
ject, provided they are massive enough to survive subsequent inter- 
actions. Given that there are no noticeable differences in the stellar 
age of the large KDCs and their host galaxies, it is likely that these 
components were formed early in these galaxies' lifetimes. Clues as 
to the formation scenario of these subcomponents may come from 
investigating their dynamical structure, using dynamical models. 



6 CONCLUSIONS 

We have presented high spatial resolution integral field spec- 
troscopy, using the OASIS spectrograph, of a subset of 28 ellipti- 
cal and lenticular galaxy centres selected from the SAURON survey. 
From these observations, we derive the distribution and kinematics 
of stars and gas, as well as stellar absorption-line indices indepen- 
dently calibrated to the Lick system. This data set has a factor four 
higher spatial sampling and almost a factor two increase in median 
spatial resolution to that of SAURON. We compare the results of 
our analysis to that of SAURON spectra within equivalent integrated 
aperture spectra. The measured parameters, including velocity dis- 
persion, emission line fluxes and stellar line-strength indices, are all 
consistent within the uncertainties, verifying our analysis. 

Within the central regions probed by OASIS, there is a rich 
variety of structure in the observed parameters, comparable to that 
observed on larger scales with SAURON. In particular, we find a 
range of stellar dynamical substructure within the central kiloparsec 
of these objects, including central disks and counter-rotating cores, 
which we quantify using the technique of kinemetry. The ionised 
gas is similarly varied, being often dynamically decoupled from 
the stellar component, and varying flux ratios indicating a range of 
ionisation mechanisms, even within the same galaxy. We analyse 
the stellar populations by comparing the observed line-strength 
maps to model predictions, giving maps of luminosity-weighted 
stellar age, metallicity and abundance ratio. These maps reveal that 
the galaxies with the youngest global populations in our sample 
also have even younger centrally concentrated populations. These 
young nuclei may be explained by young nuclear star clusters or 
other recent circum-nuclear star formation as in H II nuclei. 

The young nuclei show a propensity for harbouring compact 
kinematically decoupled stellar components, which can be associ- 
ated directly with the young stars. Considering all galaxies in the 
SAURON survey that show a clear kinematically decoupled stellar 
component, we explore the relationship between central stellar pop- 
ulation and the size of the component. We find that, although the 
'classical' decoupled components show rather old stellar popula- 
tions, there also exists a population of compact (<J 100 pc) de- 
coupled components which show a range of ages, from 15 Gyr to 
< 0.5 Gyr. Given the large difference in intrinsic size, we con- 
clude that the two types of decoupled component are not directly re- 
lated in the evolutionary sense, and that the compact variety are not 
necessarily an important feature by mass, although the young stars 
within these decoupled components are largely responsible for the 
significant spread in integrated luminosity-weighted ages observed 
for these galaxies. 

Our data suggest a link between counter-rotation and enhanced 
star-formation. By consideration of a simple model, we conclude 
that older KDCs of similar mass fractions are rather difficult to de- 
tect, and so there may be a sizeable number of low-mass counter- 
rotating components that we are missing. Likewise, co-rotating 
components of the same mass fraction leave a more subtle signa- 
ture on the observed stellar kinematics, and so can only be detected 
when the KDC population is rather young. Our sample is therefore 
too small to reliably put constraints on the importance of counter- 
rotation to recent star formation. 

These are the first observations in an ongoing follow-up of 
SAURON survey galaxies using high-spatial resolution integral field 
spectroscopy to study the central regions of these objects in detail. 
The empirical results presented in this paper are a first step towards 
understanding the intrinsic structure of the central few hundred par- 
sees of these objects, and already reveal new insights. Future appli- 
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cation of dynamical models, using both SAURON and OASIS data 
will substantiate this further. 
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APPENDIX A: MAJOR-AXIS LONG-SLIT PROFILE OF 
NGC3489 

Figure lAll presents profiles of the measured stellar kinematics, 
gas properties, line-strength indices and stellar population param- 
eters for NGC3489, extracted along a major-axis 'slit' aperture 
0'.'2 wide, collecting the value from each Voronoi bin whose centre 
falls inside. This illustrates the typical measurement error associ- 
ated with each bin in the maps presented in Fig. ll2l 



APPENDIX B: KINEMETRIC PROFILES OF STELLAR 
VELOCITY 

Fig. IB II presents profiles of kinematic position angles and peak 
absolute rotation val ues for all 28 g alaxies, measured using the 
kinemetry method of iKrainovic et alJ f2006). See Section l33l for 
details. 



APPENDIX C: COMPARISON OF OASIS AND SAURON 
SPECTRA 

Fig. lCll presents fully calibrated OASIS spectra (black lines) inte- 
grated within a 4" circular aperture centred on the galaxy nucleus 
for all the objects in our sample. Overplotted in red are equivalent 
spectra from SAURON integrated within the same 4" aperture. Both 
sets of spectra are de-redshifted, and normalized around 5100 A. 
The SAURON spectra have been broadened and resampled to match 
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Figure Al. Simulated major-axis slit aperture tf.'l wide for NGC3489. Each of the parameters derived from the two-dimensional maps are shown, with 
measurement errors, to illustrate the typical quality of our data and how the uncertainties relate to structures in the parameters. 



© RAS. MNRAS 000.17115^ 



The SAURON project VIII: OASIS/CFHT spectroscopy of SAURON E/SO galaxy centres 






t J! * 








5 1 J 1 1 
NnC3KlB 


J 








□ : 1 J J f 










> ! 





#1 ^ ] , > 

' inr^il I'll • 



Figure Bl. Radial profiles of kinematic position angle ( APA) and peak absolute rotation (ki ), measured using the kinemetry method of lKrainovic et alJ l2006h 
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the OASIS spectral resolution and sampling. Both sets of spectra 
have been 'cleaned' of gas, by subtracting an emission line model, 
as described in Section l7!4l 



APPENDIX D: COMPARISON OF OASIS AND SAURON 
EQUIVALENT WIDTH LIMITS 

Fig. lDll illustrates the difference in equivalent width limits between 
the SAURON data of Paper V and our OASIS data, for the cases 
where SAURON detects gas but OASIS does not, or where the gas 
detection with SAURON is stronger. For each galaxy, the upper pan- 
els show the estimated S/N determined from the error spectrum 
propagated through the reduction process of both data sets (circles 
= SAURON, diamonds = OASIS). This estimator is effectively an 
upper limit, as correlation effects (e.g. seeing PSF) can be difficult 
to model accurately, resulting in slightly underestimated noise esti- 
mates, at the level of ~ 5%. 

The lower panels show how these S/N values (same symbols) 
translate into equivalent width limits, using eq. (1), assuming an in- 
trinsic emission-line velocity dispersion of 50 km s~^. Overplotted 
on these limits are the measured emission-line equivalent widths 
of [OIII] from the SAURON data, where filled square symbols indi- 
cate significant detections (i.e. above the A/N threshold), and open 
square symbols are non-detections. 

This clearly shows that the poor or non-detection of emission 
with OASIS which is observed with SAURON arises simply from 
the very different 5'/A'^ of the two data sets. The lower S/N of 
the OAS I S data comes mainly from the much smaller spatial sam- 
pling on the sky, as well as generally shorter total integration times, 
smaller telescope aperture, and small difference in throughput of 
the two instruments. 

This paper has been typeset from a TgX/ KTgX file prepared by the 
author 



The SAURON project VIII: OASIS/CFHT spectroscopy of SAURON E/SO galaxy centres 57 



If TT~it~i 


IT IP ^ l^-H^ J_Ji^? 


V 11 JTGCSQS- 






^ L X ^^^^ 1 ^ 


T 




T IP J 






' LJ nOC3(iOB 


If 


■ li J1ljC4_ci_; 












L J J.I ^K^^^^ _ 

V ^ 






^ 




S.^ a.X S32 S J'i S.3t] S.:3S 




V ^ p HDC3J<i i 
9 JS B JO 9J]! SJb SJS 


S4a S.5Q S.S2 BJ4SJb BJB 



Figure CI. Comparison of OAS I S (black) and SAURON (red) spectra, plotted in In A. The spectra from both instruments have been integrated within the same 
4" circular aperture centred on the galaxy. The spectra have been de-redshifted, and normalized around 5 100 A for the purpose of this comparison. The SAURON 
spectra have been broadened and resampled to match the OAS I S system. Both data sets have been independently corrected for emission hne contamination. 
These spectra are used for the comparisons presented in Section|3] The SAURON spectra have been truncated to a wavelength range of 4825-5 280 A, to ensure a 
common wavelength coverage within the aperture for all galaxies, without encountering gaps in the field coverage due to the filter cut-off (see lKuntschner et alj 
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Figure Dl. Top panels: Compaiison of S/N as a function of circular radius for OASIS (diamonds) and SAURON (circles). Bottom panels: Corresponding 
equivalent width detection limits derived using eq.(l) for both data sets (same symbols as top panels). Overplotted with large filled (detections) and small open 
symbols (non-detections) are the [OIII] equivalent width values measured with SAURON (from Paper V). This shows that the very high S/N of the SAURON 
data results in significantly lower detection limits than that obtained with OAS I S at higher spatial sampling. This explains why some objects have clearly 
detected gas with SAURON but not with OASIS. 
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